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Abstract 

Wc examine critically the theoretical underpinnings and phenomenological implications of 
soft gluon (threshold) resummation of rapidity distributions at a hadron collider, taking Drell- 
Yan production at the Tevatron and the LHC as a reference test case. First, we show that 
in perturbativc QCD soft gluon resummation is necessary whenever the partonic (rather the 
hadronic) center-of-mass energy is close enough to threshold, and we provide tools to assess 
when resummation is relevant for a given process. Then, we compare different prescriptions for 
handling the divergent nature of the series of resummed perturbative corrections, specifically 
the minimal and Borel prescriptions. We assess the intrinsic ambiguities of resummed results, 
both due to the asymptotic nature of their perturbative expansion, and to the treatment of 
subleading terms. Turning to phenomenology, we introduce a fast and accurate method for the 
implementation of resummation with the minimal and Borel prescriptions using an expansion on 
a basis of Chebyshev polynomials. We then present results for W and Z production as well as 
both high- and low-mass dilepton pairs at the LHC, and show that soft gluon resummation effects 
arc generally comparable in size to NNLO corrections, but sometimes affected by substantial 
ambiguities. 
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1 Introduction 



Drell-Yan rapidity distributions are likely to be the standard candle which is both theoretically 
calculable and experimentally measurable with highest accuracy at hadron colliders, in particular 
the LHC The current QCD theoretical accuracy for this process is NNLO p], and small effects 
such as those related to the coupling of the gauge boson to final-state leptons have been studied 
recently [2]. Whereas the resummation of contributions related to the emission of soft gluons 
are routinely included in the computation of Drell-Yan transverse-momentum distributions 
where they are mandatory in order to stabilize the behaviour of the cross-section at low p^, their 
impact on rapidity distributions has received so far only a moderate amount of attention. This 
is partly due to the fact that even in fixed-target Drell-Yan production experiments, such as 
the Tevatron E866 jl], let alone LHC experiments, the available center-of-mass energy is much 
larger than the mass of typical final states, thereby suggesting that threshold resummation is 
not relevant. 

However, it has been pointed out since long [5] that because hadronic cross-sections are found 
convoluting a hard cross-section with a parton luminosity, the effect of resummation may be rel- 
evant even relatively far from the hadronic threshold. Indeed, in Ref. [6] threshold resummation 
has been claimed to affect significantly Drell-Yan production for E866 kinematics, though some- 
what different results have been found in Ref. [7]. It is important to observe that Drell-Yan 
data from E866 and related experiments play a crucial role in the precision determination of 
parton distributions p], so their accurate treatment is crucial for precise LHC phenomenology. 
Furthermore, threshold resummation is known [9j to affect in a non-negligible way standard 
Higgs production in gluon-gluon fusion at the LHC, even though the process is clearly very far 
from threshold. 

A systematic assessment of the relevance of soft-gluon resummation for rapidity distributions 
at hadron colliders is called for: it is the purpose of the present paper. We will take the Drell-Yan 
process as a test case, but our formalism and results can be applied to other collider processes 
such as heavy quark production or Higgs production. 

The main ingredients of this assessment are the following. First, we determine when and why 
threshold resummation is relevant. Then, we discuss how the resummed contribution is defined, 
and specifically we deal with the divergence of the series of resummed terms. Next, we address 
the issue of combining resummed and fixed-order results. Finally, we turn to phenomenology 
and assess the size and theoretical uncertainty of resummation at the Tevatron and LHC. Each 
of these steps turns out to be nontrivial, as we now briefly sketch. 

The standard physical argument to explain why resummation may be relevant even when 
the hadronic process is relatively far from threshold goes as follows [lOJ. The quantity which is 
resummed in perturbative QCD is the hard partonic cross-section, which depends on the partonic 
center-of-mass energy and the dimensionless ratio of the latter to the final state invariant mass. 
Therefore, resummation is relevant when it is the partonic subprocess that is close to threshold. 
The partonic center-of-mass energy in turn can take any value from threshold up to the hadronic 
center-of-mass energy, and its mean value is determined by the shape of the PDFs: therefore, 
one expects threshold resummation to be more important if the average partonic center-of-mass 
energy is small, i.e. if the relevant PDFs are peaked at small x (such as gluons or sea quarks, as 
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opposed to valence quarks). This for instance explains why threshold resummation is especially 
relevant for Higgs production in gluon-gluon fusion. 

We will show that this can be made quantitative using a saddle-point argument in Mellin 
space: for any given value of the hadronic dimensionless variable r the dominant contribution 
to the cross section comes from a narrow range of the variable N which is conjugate to r upon 
Mellin transform. In Mellin space the cross-section is the product of parton distributions (PDFs) 
and a hard coefficient, but it turns out that the position of the saddle is mostly determined by 
the PDFs. Moreover, the result is quite insensitive to the non-perturbative (low-scale) shape 
of the parton distribution and mostly determined by its scale dependence, specifically by the 
low-x (or low-A^) behaviour of the relevant Altarelli-Parisi splitting functions: the faster the 
small-x growth of the splitting function, the smaller the average partonic center-of-mass energy, 
the farther from the hadronic threshold the resummation is relevant. This is reassuring, because 
it means that the region of applicability of threshold resummation is controlled by perturbative 
physics. The issue of the persistence of sizable soft gluon emission terms even far from threshold 
was also addressed, using methods of soft-coUinear effective theory, in Ref. [7j, where in the large 
T > 0.2 region it was related to the (non-perturbative) shape of parton distributions, though it 
was also observed for smaller r values. Q 

Having established the region in which threshold resummation is relevant, we have to face the 
fact that resummation sums an infinite series of contributions to the expansion of the hard par- 
tonic cross-section in powers of the perturbative strong coupling ag(Q^) (with the hard scale 
of the process) which diverges if the resummation is performed at any given logarithmic order. 
This divergence can be treated in various ways: here we will consider the minimal prescription 
(MP) [12] and the Borel prescription (BP) [131 E]) in both of which the resummed result is 
a function of as{Q^) to which the perturbative expansion in powers of as{Q^) is asymptotic. 
Both the BP and MP can be obtained by adding to the perturbative series a contribution which 
removes the divergence, which is power-suppressed in for the BP, and it has support in an 
unphysical kinematic region (below threshold) for the MP. 

In practice, however, different resummation prescriptions differ not only because of the way 
the high-order divergence of the expansion is handled, but also because if they are applied to a 
low-order truncation of the divergent series they differ by subleading terms. We will show that 
in fact this difference turns out to be by far phenomenologically the most significant, unless one 
is very close to the hadronic threshold, which is in practice a very rare occurrence. We will thus 
assess the ambiguity that subleading terms induce on resummed results, and to which extent 
they can be optimized in order to ensure stability of the resummed expansion once resummed 
results are matched to the standard fixed-order ones. 

With one (or more), possibly optimal, resummation prescriptions at hand, we turn to the 
resummation of rapidity distributions. It turns out that this can be performed by relating 
the resummed expressions to those of the inclusive cross-section. This entails some further 
ambiguities in the treatment of subleading terms, though we shall see that these are in practice 
phenomenologically very small. Finally, we will implement the resummation up to the next-to- 

"'^ Shortly after this paper appeared in preprint form, the treatment of this issue in soft-collinear effective theory 
was revisited in a quantitative way in Ref. where it was related to a parameter determined by the shape of 
parton distributions. 
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next-to-leading log (NNLL) level combined with the next-to-next-to-leading fixed-order (NNLO) 
result. We will show that this can be done efficiently by projecting resummed results on a basis 
of Chebyshev polynomials: with this approach, resummed results can be easily obtained using 
any external set of parton distributions. This will enable us to obtain predictions at the Tevatron 
and LHC. We will see that even at the LHC the impact of the resummation is not negligible, 
and comparable to the size of the corresponding fixed-order corrections, especially in the central 
rapidity region. 

The structure of this paper is the following: in Section 2 we will present in detail the ar- 
guments which allow one to determine the relevant partonic center-of-mass energy for given 
hadronic kinematics and parton distribution, and thus to assess the relevance of threshold re- 
summation. In Section 3 we will discuss and compare the minimal and Borel prescription for 
resummation, the way they can be matched to fixed-order expressions, the way subleading terms 
are treated with the various prescriptions, and the associate ambiguities. In Section 4 we will 
present the general formalism for the construction of resummed expressions for rapidity dis- 
tributions from their inclusive counterparts, and the numerical implementation of resummed 
results. In Section 5 we will turn to phenomenology for the production of neutral and charged 
Drell-Yan pairs, with different values of their invariant mass at Tevatron and LHC energies. We 
will compare theoretical predictions both to E866 and recent CDF data. 

2 When is threshold resummation relevant? 

As we have mentioned in the introduction, both simple physical arguments and evidence from 
explicit calculations [10] suggest that threshold resummation may be relevant even quite far 
from the hadronic threshold, provided the partonic average center-of-mass energy is sizably 
smaller than the available hadronic center-of-mass energy. Here we derive this conclusion from 
a quantitative argument, which will allow us to assess the relevance of threshold resummation 
for a given process. We assess the impact of parton distributions by means of a Mellin-space 
argument. For a given process and a given value of r = Q'^/s, we determine the region of 
the variable N, Mellin-conjugate to r, which provides the dominant contribution to the cross- 
section. First, we show that the saddle point is mostly determined by the small-x behaviour of 
the PDFs, which in turn is driven by perturbative evolution. Next, we determine the impact of 
the inclusion of PDFs for the Drell-Yan process. Finally, with specific reference to the Drell-Yan 
process, we assess the N region where threshold resummation is relevant. 

2.1 The impact of PDFs: A^— space vs. r— space 

The cross-section for a hadronic process with scale and center-of-mass energy s = Q'^/t can 
be written as a sum of contributions of the form 




(2.1) 
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in terms of a partonic cross-section a and a parton luminosity, in turn determined in terms of 
parton distributions fi{ 



C{z)= C—h{x,)h[^]. (2.2) 



Here we denote generically by a a suitable quantity (in general, process-dependent) which has 



the property of factorizing as in Eq. (2.1). Such quantities are usually related in a simple way 



to cross-sections or distributions; for example, in the case of the invariant mass distribution of 
Drell-Yan pairs, a is in fact ^^i- 



In general, the cross-section gets a contribution like Eq. (2.1) from all parton channels which 
contribute to the given process at the given order, but this is inessential for our argument, so 
we concentrate on one such contribution. 



In Eq. (|2.1|), the partonic cross-section, which is computed in perturbation theory, is evaluated 

r 

X1X2S, (2.3) 



as a function of the partonic center-of-mass energy 



T / Z 

where X\ and X2 = zjx\ are the momentum fractions of the two partons. Therefore, the threshold 
region, where resummation is relevant, is the region in which s is not much larger than . 
However, all values of X\^X2 between r and 1 are accessible, so whether or not resummation is 
relevant depends on which region gives the dominant contribution to the convolution integrals 



Eqs. (2.1, 2.2). This dominant region can be determined using a Mellin-space argument. 
The Mellin transform of (T(r, Q^) is 

a(iV,Q2)= [\Tr''-'a{T,Q'), (2.4) 

with inverse 

1 pN+ica -| pN+ica 

a(r,Q2) = ' / dNT-''a{N,Q') = ^ / rfiV e^(^'^^«'), (2.5) 

27rZ Jtv-ioo JN-ioo 

where N is larger than the real part of the rightmost singularity of (t{N, Q"^) (by slight abuse of 
notation we denote with a both the function and its transform), and in the last step we have 
defined 

E{t,N;Q^) = N\n- + \na{N,Q^). (2.6) 
r 

The function (t{N, Q"^) has a singularity on the real positive axis because of the parton luminosity; 
to the right of this singularity, it is a decreasing function of N , because the area below the curve 
r^~V(r, Q^) obviously decreases as N increases. As a consequence, E{t, N; Q"^) always has a 
minimum on the real positive axis at some = A^o- Hence, the inversion integral is dominated 
by the region of A^ around Nq, and can be approximated by saddle-point, expanding E{t, N; Q"^) 
around A'o- 
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When r — )■ 1, the slope of the straight hne A^ln^ decreases, and the position A^o of the 
minimum is pushed to larger values, so the large-r behaviour of cr(r, Q"^) is determined by the 
large- behaviour of a{N,Q'^), as it easy to show using the saddle point approximation: in 
fact, for typical parton distributions in this limit the saddle point approximation becomes exact. 
This is shown in Appendix |A| where some properties of the saddle point approximation to Mellin 
transforms are collected. 

The position of the saddle point A'o is strongly influenced by the rate of decrease of the cross- 



section (t(A^, Q^) as A^ grows. Indeed, in Mellin space, the cross-section Eq. (2.1) factorizes: 



a{N, Q') = CiN, Q') aiN, asiQ')). (2.7) 

It is then easy to see that the decrease of (j{N, Q"^) with A^ is driven by the parton luminosity 
jC,{N,Q'^): in fact, it turns out that, for large A^, a{N,as{Q'^)) is an increasing function of A^. 
Indeed, it turns out (see Appendix |A]) that a distribution, unlike an ordinary function, need not 
be a decreasing function of A^. However, the parton luminosity always offsets this increase if 
the convolution integral exists, because the cross-section a{T, Q"^) is an ordinary function. For 
example, the 0{as) Drell-Yan partonic cross-section is shown in Fig. [l| where it is clear that 
even though at small A^ the cross-section decreases with A^, at large A^ it increases. 



Drell-Yan partonic q-qbar. Order Mellin transform 




Figure 1: The 0{as) neutral current Drell-Yan partonic coefficient Ci{N), Eq. (2^25), plotted as 
a function of A^ (solid curve). The logarithmic approximations C|°^(A^) Eq. ( |2.27 ) and C}°^ (A^) 
Eq. (2.29) are also shown (dotted and dashed curves, respectively). 



As a consequence, when r is large, the position of the saddle point A^o is completely controlled 
by the drop of the parton luminosity: indeed, in the absence of parton luminosity, the saddle 
point would be very close to the minimum of the partonic cross-section, which is around A^ ^ 1. 
When T is smaller, even without PDF the location of the saddle is controlled by the partonic 
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cross-section, which in this region is a decreasing function of A^. However, this decrease is 
much stronger in the presence of a luminosity, so the location of the saddle is substantially 
larger. Hence, in the large r region the effect of the resummation is made much stronger by the 
luminosity, while for medium-small r if the luminosity decreases fast enough, Nq may be quite 
large even if r ^ 1, i.e. far from the hadronic threshold, thereby extending the region in which 
resummation is relevant. 

The position of the saddle point in the various regions can be simply estimated on the basis 
of general considerations. At the leading log level, parton densities can be written as linear 
combinations of terms of the form 



/,(iV,Q2)=exp 



7.(iV) j^ «s(gg) 



/3o '"«s(Q2) 

where 7j(A^) are eigenvalues of the leading-order anomalous dimension matrix: 



din f,{N,Q^ 
dlnQ^ 



as{Q'H{N) + Oialy, 



(91nQ2 



(2.9) 



and f-^\N) = fi{N,Ql) are initial conditions given at some reference scale Q^. The cross- 
section is correspondingly decomposed into a sum of contributions, each of which has the form 



of Eq. (2.5), with 



+\nf\'\N) + \nff\N) + \na{N,as{Q^)). 



E{t,N;Q' 

(2.10) 

At large A^, this expression is dominated by the first term, which grows linearly with A^, while 
at small N the behaviour of E{t, N; Q"^) is determined by the singularities of the anomalous di- 
mensions, which are stronger than those of the initial conditions if > Qq, given that low-scale 
physics is both expected theoretically from Regge theory [151 US] and known phenomenologically 
from PDF fits [IT] to produce at most poles but not essential singularities such as those obtained 
exponentiating the anomalous dimensions. Indeed, assuming a power behaviour for f-^^ (z) both 
at small and large z, 



^(l-z) 



SO that 



In Z^'''''' ( A^) behaves as InA^ both at large and small A^, and is thus subdominant in comparison to 



fr\N) 



r(Ar + a,)r(A + i) 

T{N + a, + A + 1) 



(2.11) 



(2.12) 



either the r dependent term or the anomalous dimension contribution in Eq. (2.10). A similar 
argument holds for the partonic cross-section term lno"(A^). 

The position of the minimum is therefore determined by the transition from the leading 
small-A^ drop due to the anomalous dimension term and the leading large-A^ rise due to the 
r-dependent term, up to a correction due to the other contributions to Eq. (2.10). When r is 



large, the rise in the first term is slow, and it only sets in for rather large A^ so the correction 
due to the other contributions may be substantial. This is the region in which resummation is 
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surely relevant because the hadronic r is large. But when r is not so large, the rise sets in more 
rapidly, in the region where the second term is dominant and the correction from the initial 
PDFs and the partonic cross-section is negligible. 

In order to provide an estimate of the value of A^o in either region, we note that the leading- 
order Altarelli-Parisi anomalous dimensions at small N behave as 



1 



7i 



+ less singular. 



(2.13) 



where "less singular" denotes terms whose singularity has a smaller real part, and A'p = for 
Iqqilqgi and Np = 1 for ■jgqy'jgg- Thls implies that one of the two singlet anomalous dimension 
eigenvectors has a leading singularity at = 1, while the other has a singularity at = and 
it is thus suppressed at small A^, like the nonsinglet anomalous dimension which is thus also 
suppressed. This pattern persists to all perturbative orders. It follows that singlet quark and 
gluon distributions have a steeper small-A^ and thus small-z behaviour. 

Expanding the anomalous dimension about its rightmost singularity at leading order we have 



7+(iV) 



(0) 

7+ 
A^- 1 



[1 + 0{N-1)] 



7ns(iV) 



/ns 



[1 + om 



(2.14) 



where 7+ and 7ns are respectively the dominant small-A^ singlet eigenvalue and nonsinglet 
anomalous dimension, and 

- (2.15) 



(0) 

7V 



^(0) 

/ ns 



(because of momentum conservation), while 7ns(l) = (because of 



We note that 7+ (2) 

baryon number conservation): hence, we expect the small A^ approximation Eqs. (2.14, 2.15) to 
break down around N ^ 2 — ki, with 



0; 



1. 



(2.16) 



We can then consider three cases, according to whether 7j, 7^- in Eq. (2.10) are both singlet, 
both nonsinglet, or one singlet and one nonsinglet. The three cases correspond respectively to 
the leading behaviour of, for instance, Higgs production in gluon fusion, Drell-Yan production 



at the Tevatron, and Drell-Yan production at the LHC Substituting Eq. (2.14) in Eq. (2.10) 
and neglecting the last three terms the saddle is found to be at 



1 kikj -\- 



\/3olni as(Q2 



(2.17) 



where i, j take the values ns and +, ki are given by Eq. (2.16), and 



7 



(0) 



2^(0). 

'ns ' 



7?! = 27f ; 



(0) 
7+ns 



(0) 

7+ • 



(2.18) 



We expect Eq. (2.17) to provide a good approximation to the position of the saddle point in the 
region 1 — kikj < N <^ 2 — kikj, while for larger values of A^ the last three terms on the right-hand 
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side of Eq. (2.10) provide a correction of increasing size. This correction cannot be estimated in 
a universal way, in that it will generally depend on the shape of the initial parton distributions 

f^^\ However, note that at large both the nonsinglet and singlet anomalous dimension drop: 

r(jv) 



they are negative, and their modulus grows as InA^. Furthermore, = A^-^(l + 0{1/N)), 

which, together with Eqs. (2.11, 2.12) implies that also the contributions from ff^^ to Eq. (2.10) 



drop logarithmically as — ?■ oo. This drop is a generic feature of the initial PDFs: it is 
necessary in order for the convolution integral which gives the physical cross-section to exist (see 
Appendix |A]); likewise, the drop of the anomalous dimension (at most as a power of InA^) is a 
generic feature of perturbative evolution, related to momentum conservation of gluon radiation, 
and it persists to all orders. 



It follows that at large the 1/A^ drop Eq. (2.14) is replaced by a slower logarithmic drop, 
with a coefficient determined by both the large-iV behaviour of the anomalous dimension and 
by the initial PDFs. This leads to a rather larger value of the saddle point Nq than given by 
Eq. (2.17). The value of A^°- as a function of r given by Eq. (2.17) is plotted in Fig. |2]for the 
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Figure 2: The position of the saddle point Nq for the Mellin inversion integral Eq. (2.5) with 



the exponent Eq. (2.10) as a function of r determined neglecting a and using toy initial PDFs 
Eqs. (2.11, 2.19) and leading-order anomalous dimensions. Upper curves: exact value; lower 
curves: the approximation Eq. (2.17). The curve shown is obtained with Qo = 1 GeV, Q = 
100 GeV. 



three cases of Eq. (2.18): in each case, we show both the result obtained using the approximation 
Eq. (2.17), and that found by exact numerical minimization of Eq. (2.10), with the full leading- 
order expression of the pertinent anomalous dimensions 7i, 7j, 



and initial PDFs of the form 
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Eq.([2lTl), with 



Or, 



/3n 



0; (3, 



(2.19) 



but still neglecting the \na{N,as{Q^))- We see that indeed the pole approximation Eq. (2.17) 
holds reasonably well for 1 — kikj < N <^ 2 — kikj. 

We can get a feeling for the implication of this for the r region in which resummation is 
relevant by very roughly taking the value = 2 as that of the transition from the small-A^ 
region to the large- region. Indeed, this value corresponds to the energy-momentum operator 
and it thus sets the transition between the large N region in which parton distributions drop 
(while growing in modulus) and the small N region in which they grow as the scale is raised. 
This rough estimate is borne out by a more quantitative assessment of the transition point which 



will be performed in Sect. 2^ We then see that, if at least one of the parton distributions is 
flavour singlet, the position of the saddle point remains in the large- region down to fairly low 
values of r ~ 0.01. This extension of the region where resummation is relevant to small r is due 
to the rise of the anomalous dimension related to the pole at = 1 Eq. (2.13) in the singlet 
sector. 

Larger values of r ~ 0.1 correspond to values of Aq which are anyway in the large- A^ region. 
However, in this region the rapid drop of the parton distribution, due both to its initial shape 
and to its evolution, greatly increases the impact of resummation by raising the position of the 
saddle, which at r ~ 0.2 is already A'o ~ 6. 

2.2 The impact of PDFs: the Drell-Yan process 

We now assess the impact of parton distributions in the specific case of Drell-Yan production. 



In this case, the quantity cr(r, Q"^) which appears in Eq. (2.1) is given by 



-ir,Q^) = ~ir,Q% (2.20) 

where ^^5^(t, Q^) is the invariant mass distribution of Drell-Yan pairs. The corresponding 
partonic quantity cr{z,as{Q^)) of Eq. (2.1) is then 



a{z,as{Q')) 



1 daoY 



:^,«s(Q')). 



(2.21) 



z dQ^ 

It is further convenient to define a dimensionless coefficient function C{z,aa,{Q'^)) through 

a(z,as(Q')) = <^oC(z,as(Q')). (2.22) 
The coefficient function admits the perturbative expansion 



C(z,as) 



5{l-z) + '^C^{z)+(^)' C2{z) + 

TT V TT / 



(2.23) 



the Born term (Tq contains the electroweak coupling and it is given e.g. in Ref. [18]. Threshold 
logarithms only appear in the quark-antiquark channel, while in other partonic channels they 
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Figure 3: The position of the saddle point A^o for the Melhn inversion integral Eq. (2.5) as a 



function of r with the cross-section Eq. (2.7) determined using the 0{as) Drell-Yan cross-section 
Eq. (|2.25|) for neutral dileptons and NNPDF2.0 g] parton distributions, with Q = 100 GeV. 



The two upper curves refer, from the top, to pp and pp collisions. The lowest (dashed) curve is 



the position of the saddle at the parton level, i.e. omitting the parton luminosity in Eq. (2.7) 



are suppressed by powers of 1 — z. In the quark-antiquark channel the next-to-leading order 
contribution is given by [19] 



ln(l - z) 



4 1 — z 

In v^- 2(1 + z) In 



1-z 



71 

y 



4)6{l-z)\. (2.24) 



Note that of course at leading order cr{N) is just a constant. The Mellin transform of the NLO 
term is 



4 + 27^ + 2^1{N) - MN) +MN + 2)+ Aj^MN) 



+ ^[7. + ^o(iV + l)] + ^ 



[^^ + MN + 2)] }, (2.25) 



and was shown in Fig. [TJ 

We have then determined the position of the saddle point Nq in a realistic situation, i.e.. 



using the partonic cross-section Eq. (2.25) for Drell-Yan production of a neutral lepton pair 



of invariant mass Q = 100 GeV at a pp or pp collider, with a parton luminosity determined 
using NNPDF2.0 [8] parton distributions. Comparing the realistic curves of Fig. [3] to those in 
Fig. [2| which were determined using the toy initial PDFs Eqs. (2.11, 2.19) and neglecting the 
contribution from the hard cross-section we see that the pp curve of Fig. |3] agrees well with the 
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sea-valence curve of Fig. |2j as one would expect given that in pp collisions one must always 
pick up at least a sea (antiquark) PDF. The case of pp is slightly more subtle: in this case, for 
T > 0.1 the curve in Fig. [3] agrees well with the valence-valence curve of Fig. [2] However, for 
smaller values of r the position of Nq computed using the full luminosity decreases much more 
slowly: this is due to the fact that as < 2 the contribution 7+ rapidly grows due to the pole 
Eq. (2.14) so that even the valence distribution is dominated by it. As a consequence, even in 
the pp the resummation region is further extended to somewhat lower r values than it would be 
the case for a pure valence-valence luminosity. 

In Fig. |3] we also show the position of the saddle that is obtained if the parton luminosity 
is omitted, i.e. for the parton-level cross-section. It is clear that, as we argued in Sect. |2l2, the 
position of the saddle is determined by the PDFs, and it is much larger than that found at the 
parton level, thereby supporting the conclusion that the convolution with the parton luminosity 
greatly enhances the importance of resummation, and it extends it to a much wider kinematic 
region. 

In summary, we conclude that the resummation region N > 2 corresponds to r > 0.003 for 
pp collisions, and r > 0.02 for pp collisions. For r < 0.1 the position of the saddle is determined 
by the pole in the anomalous dimension, while for larger values of r the large x drop of PDFs, 
due both to their initial shape and to perturbative evolution, very substantially enhances the 
impact of resummation. These values have been obtained for \/Q^ = 100 GeV. We have checked 
that they depend very weakly on Q^, which is expected, because enters the determination of 
the saddle point only through the scale dependence of the parton densities. 



2.3 The resummation region for the Drell-Yan process 



Having established that the PDFs drive the position of the saddle point, and having explicitly 
determined this value, we now have to determine quantitatively the value of N at which resum- 
mation becomes important. So far, we have taken conventionally the momentum-conservation 
point iV = 2 as the value at which logarithmically enhanced contributions give a sizable contri- 
bution to the cross-section. We would now like to establish this in a quantitative way. 

To this purpose, we compare Ci{N) Eq. (2.25) to its logarithmic approximation. The 
logarithmically-enhanced term in Ci{z) Eq. (2.25) is 



ln(l 



whose Mellin transform is 



TT 



+ 27' 



(2.26) 



(2.27) 



The function Ci (N) was displayed as a function of along the real positive axis in Fig. [T] The 
logarithmic approximation Eq. (2.27) is also shown in the same figure. It is clear that the full 



result and the log approximation agree at large A^, up to a small constant shift. Indeed, 



lim 



Ci{N)-C{°^{N) 



TT 



-4 



(2.28) 
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For N slightly above 2 the logarithmic contribution is already about 50% of the full result, but it 
rapidly deviates from it as decreases. This suggests that indeed the logarithmic contribution 
is sizable for N '>2. 

It is interesting to observe that there is a certain arbitrariness in the definition of the log- 



arithmically enhanced contribution: indeed, the Mellin transform Eq. (2.27) of the logarithmic 



contribution Eq. (2.26) contains terms which are not logarithmically enhanced in iV-space: the 
Mellin transform of the next '^-leading-log x expression only coincides with the next*^-leading-log 
expression up to subleading terms. For example, ipi ~ ^, so that we could equally well 



identify the logarithmically enhanced contribution to Eq. (|2.25|) with 



2^/>o^(iV)+47,Vo(iV) 



TT 

y 



+ 2j 



which is the Mellin transform of 



{ 


"ln(l - z)' 


Inv^j 




_ I- z _ 


+ 1-4 



(2.29) 



(2.30) 



In other words, logarithmically enhanced contributions in A^-space also contain subleading terms 
when transformed to z-space, and conversely (see Appendix [B| for details on the relevant Mellin 
transforms) . 

The curve corresponding to C}°^ (A^) is also shown in Fig. [l| It is clear that in the large 
N >2 region it differs from C^°^{N) by a negligible amount. However, at small A^ this particular 
logarithmic approximation turns out to be much closer to the full result. Whereas of course the 
choice of subleading terms to be included in the resummation is essentially arbitrary, in practice 



some choices may be less natural than others: this will be discussed in Sect. |3.2| below. 

In order to check that these conclusions are generic, we repeat the comparison at the next 
perturbative order. The order Drell-Yan cross-section 6*2(2) has been computed in Ref. pO] . 
In Fig. [4] it is compared to its logarithmic approximation C^2^{N), defined as the Mellin transform 
of all contributions to C2{z) which grow logarithmically as z — )■ 1, as well as another logarithmic 
approximation C^2^ (N) which differs by the former through power-suppressed terms, analogous 
to C|°^ (A^) and to be discussed in more detail in Sect. 



3.2 



below. 



It is apparent from Fig. |4] that the conclusions drawn at NLO are qualitatively unchanged, 
though there is some quantitative difference: the logarithmic approximation, which differs by a 
constant from the asymptotic large- A^ behaviour of the coefficient function, provides a sizable 
or even dominant contribution to it for A^ > 2 — 3, according to the choice of subleading terms: 
approximations which differ by subleading terms are all very close for A^ ^ 3, though they may 
differ substantially at small A^. 

As a final comment we note that the region in which we have found logarithmic effects to 
be important is in fact rather wider than the region in which as In A^ ~ 1. In this region even 
though logarithmically enhanced terms may lead to a substantial contribution, they behave in 
an essentially perturbative way, in that (as In^ A^)'^^^ < (asln^ A^)'^, and therefore the all-order 
behaviour of the resummation is irrelevant. In this intermediate region, the resummation may 
have a significant impact, but with significant ambiguities related to subleading terms. 
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3 Resummation 

After having provided an assessment of the kinematic region in which logarithmically enhanced 
terms are significant, we now discuss their all-order resummation. We will assume knowledge of 
the resummed results in A^-Mellin space, both at the level of inclusive cross-sections [2T1 122| [23] 
and rapidity distributions [211 ESI [6], and we will concentrate on the ambiguities related to 
their definition. First, we briefly review some prescriptions introduced in order to deal with the 
divergent nature of the perturbative expansion for resummed quantities. Next, we discuss the 
role and impact of subleading terms in the resummation, with specific reference to the subleading 
terms which are introduced by different resummation prescriptions. 



3.1 Resummation prescriptions 



Resummation of the partonic cross-section ci{z, as{Q )) Eq. (2.1 ) is most naturally performed in 
space, because the solution of the relevant Altarelli-Parisi equation [22j, the relevant factor- 
ization theorem [2T1 126] and the relevant renormalization group equation [26l [23] are all naturally 
formulated in A^-space. The underlying physical reason is that, in the soft limit, while ampli- 
tudes factorize due to the eikonal approximation both in A^- and in z-space, the relevant phase 
space only factorizes in A^-space but not in 2;-space [27] . 

Recently, resummation has also been performed using SCET techniques both in A^-space [28] 
and in momentum space [291 |7j; in the latter approach the soft scale whose logarithms are re- 
summed is not the partonic scale Q^{1 — z) but rather a soft scale fis, independent of the parton 
momentum fraction z, but related to the hadronic scale Q^{1 — t). As a consequence, the hard 
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coefficient function depends on r not only through the convolution variable, but also directly 
through the soft scales: therefore, the resummed result can no longer be factorized by Mellin 
transformation into the product of a parton density and a hard coefficient. For this reason, a 
direct comparison between the result obtained through the SCET approach and that based on 



standard factorization Eq. (2.1) is not possible at the parton level. A phenomenological com- 



parison is possible [7j, but it requires either assuming a specific form of the parton distribution 
functions, or switching to a SCET formulation in which fig depends on A^, which however is not 
the Mellin transform of the r-space SCET result, advocated for phenomenolgy in Ref. [7]. 

The general structure of A^-space resummed expressions can be understood by considering 
the case of the inclusive coefficient function C{N,as), obtained by Mellin transformation from 
C{z,as) Eq. (2.23). At the resummed level. 



a 



■ giW + g2{X) + a gsW + ^4(A) + 



(3.1) 

(3.2) 
(3.3) 



where a is a pro cess- dependent constant (a = 2 for Drell-Yan production). At the next'^- 
to-leading logarithmic (N'^LL) level functions up to gk+i must be included, and go must be 
computed up to order (see Appendix |d] for explicit expressions). Expansion of the resummed 
coefficient function Eq. (3.1 ) in powers of as(Q^) up to order n gives the logarithmically enhanced 



contributions to the fixed-order coefficient functions Ci{N) Eq. (2.23) up to the same order, 
with, at the N*^LL level, all terms of order In™ with 2{i — k) < m < 2i correctly predicted. 
The inclusion up to the relevant order of the function go is necessary, despite the fact that go 
does not depend on InA^, because of its interference with the expansion of the exponentiated 
logarithmically enhanced functions gi with i > 1. 

Matched resummed coefficient functions are then obtained by combining the resummed result 
Eq. (3.1) with the fixed-order expansion in power of as, and subtracting double-counting terms, 
i.e. the expansion of C^'^^{N,as{Q'^)) in powers of as{Q^) up to the same order: 



^N'^LL+NPLO 



(AT, as) = C-{N, as) + E (^)' C^iN) - 

i=o j=o 



ai 



d^C'^'iN, as) 
dai 



(3.4) 



as=0 



However the perturbative expansion of C'^'^^{z,as{Q'^)) Eq. (3.1) in powers of as(Q^) turns 
out to be divergent. This follows from the fact that the functions gt in Eq. (3.2) depend on A^ 
through 



c^siQVN^) 



as{Q' 



{1 + 0{as{Q'))) ; L = a\n 



(3.5) 



with a defined 
cross-section in 

a{N,asiQ^)) at 
axis from Nl = 
C^^^iN,asiQ')) 



1 + L 

in Eq. (3.3). As a consequence, the expansion of the resummed partonic 
powers of asiQ"^) has a finite radius of convergence dictated by \L \ < 1, and 

fixed asiQ"^) has a branch cut in the complex A^-plane along the positive real 
exp(l/a;) to +00. But a Mellin transform always has a convergence abscissa, so 
cannot be a Mellin transform. 
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On the other hand, any finite-order truncation of the series expansion 



C^^%N,as{Q')) = J2 



TT 



(3.6) 



behaves as a power of IniV at large and it is thus free of singularities for large enough. 
Hence, C^^^{N, ag(Q^)) can be viewed as the Mellin transform of the function C^^^{z, asiQ"^)) 



TT 



(3.7) 



such that 



N+ioo 



N- 



dN 
2Tii 



-N 



^res 



(iV). 



(3.^ 



It follows, by contradiction, that the series Eq. (3.7) must diverge fT3]. It turns out [13] that 



if the Mellin inversion Eq. (3.8) is performed to finite logarithmic accuracy, the series Eq. (3.7) 



acquires a finite but nonzero radius of convergence in z] however, this does not help given that 



the convolution integral Eq. (2.1) always goes over the region where the series diverges. 



Hence, any resummed definition must either explicitly or implicitly deal with the divergence 



of the perturbative expansion Eq. (3.7). We now consider two prescriptions in which this is done 



by constructing a resummed expression to which the divergent series is asymptotic. 



3.1.1 Minimal prescription 

The minimal prescription (MP) [12] defines the resummed hadronic cross-section as 



o'mp(t, 



c+ioo 



dN T-''C{N, Q') ^■•''^(iV, as(Q')), 



(3.9) 



The integration path in Eq. (3.9) is taken to the left of the 

that the cross-section 



where (Jq is defined in Eq. (2.22) 

cut, but to the right of all other singularities. It is shown in Ref. [12[ 
obtained in this way is finite, and that it is an asymptotic sum of the divergent series obtained 



by substituting the expansion Eq. (3.6) in Eq. (3.9) and performing the Mellin inversion order 
by order in as(Q^). Of course, if the expansion is truncated to any finite order the MP simply 



gives the exact inverse Mellin transform, namely, the truncation of Eq. (3.7) to the same finite 
order. 

If the MP is applied to the resummed partonic cross-section, i.e. if one omits the parton 
luminosity C{N,Q^) in Eq. (3.9) then, because of the branch cut to the right of the integration 
contour, the ensuing integral gives a function (Jmp(-z, Q^), which does not vanish when z > 1 [12]. 
However, as shown in Ref. \12l, the contribution from the z > 1 region is exponentially suppressed 



in . In the vicinity of 2; = 1 the integral oscillates strongly, as shown in Fig. 
o'Mp{z, Q"^) is displayed for the Drell-Yan resummed cross-section, evaluated at two seal 
are relevant for the phenomenological discussion of Sect. |5} When folded with the luminosity, the 
hadronic cross-section receives a contribution from the unphysical region and the corresponding 



where 
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Figure 5: The partonic cross-section (TMpfQ ^, z) computed using the minimal prescription, 



i.e. setting C{N,Q ) = 1 in the integral Eq. (3.9h, evaluated at 



8 GeV and 



100 GeV. The curve shown is obtained using on the r.h.s. of Eq. (3.9) the NLL expression for 



the Drell-Yan coefficient function Eq. (3.1, 3.2) with only gi and g2 included. 



integral thus does not have the form of a convolution. The ensuing integral is finite, but the 
oscillatory behaviour of the partonic cross-section makes its numerical computation difficult. A 
technical solution to this problem is provided in Ref. |12]; here, we will propose in Sect. |4] a 
different solution, and use it for phenomenology in Sect. [5j 

3.1.2 Borel prescription 

An alternative prescription is based on the Borel summation of the divergent series. This pre- 
scription was developed in Refs. [131 E]; here we give an equivalent, but a somewhat simpler 
presentation of it. To this purpose, it is convenient to rewrite the resummed coefficient function 



C-^(Ar,as(Q2)) Eq. (3.1) as 



C^'^N, a,{Q^)) = 1 + a,{Q')), (3.10) 
where Cq^^{N, ctg) = 1 is the (iV-independent) Born result, so that only logarithmically enhanced 



terms are included in E. Using Eq. (|B.6|) we see that 

S(^,«3(Q')) 



N+ioo 



2ni jN^ioo 

oo 

R{z) = Y,hk{as{Q'))a'ck{z), 



Ini 



k=l 



where 



In^ 



and /ife(Q;s(Q^)) are the coefficients of the expansion 

oo 



(3.11) 



(3.12) 



(3.13) 



k=l 
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whose as((5^) dependence will henceforth be omitted for notational simplicity. Using Eq. (B.5) 
Eq. ( 3.11[ ) can be written as 



Riz) 



1 



In^ \ 



k=l 



(3.14) 



The Borel prescription can be now formulated. First, we note that the Borel transform of 
R{z) with respect to a, which is found replacing 



w 



{k-l)\ 

in Eq. (3.14), is convergent, and can be summed in closed form: 



R{w, z) 



1 

27ri 



ln« 1 d 



i r(0 dw \i 



w 



,«s(Q' 



(3.15) 



(3.16) 



fc=i 

The branch cut of S(L,as(Q^)), — oo < L < —1, is mapped onto the range —w < ,^ < on 
the real axis of the complex ^ plane. Hence, the ^ integration path is any closed curve which 
encircles the cut. 

Next, we observe that the inverse Borel transform oi R[w, z) does not exist because it involves 
evaluation of the function l/r(^) on the negative real axis where it is badly behaved as ^ — )■ 
— oo [13]. We may however define the inversion integral by introducing a cutoff at some finite 
value C of the Borel variable w: if we replace R{z) Eq. (3.14) with 



c 



dw e "R{w,z) 



1 

27d 



di ln« \ 



e r(0 

then the resummed function S(2;, a;s(Q^)) becomes 



c 



dw e 



— : 

dw 



w 



«s(Q') 



Ec(z,«s(g')) 

It is proved in Refs 



1 

27x1 



d^ 1 



c 



dw e 



— : 

dw 



,«s(Q' 



ln« 



(3.17) 



(3.18) 



that the original divergent series Eq. (3.11 ) for S is asymptotic to the 
function Sc(z, as(Q^))) and furthermore, that for any finite-order truncation of the divergent 
series, the full and cutoff results differ by a twist- (2 + ^) term, where C is the cutoff. The 
parameter C can be chosen freely in the range C > a (with a as in Eq. ( 3.3[ )), with different 
choices differing by power suppressed terms. We will use henceforth the "minimal" choice C = a. 

A somewhat simpler result is obtained if the Borel transform is performed through the 
replacement 



1 w^ 
a k\ 



in Eq. (3.14), instead of Eq. (3.15). In this case, instead of Eqs. (3.17, 3.18) one gets 

di ln« \ 



Rc{z) 



Sc(^,«s(Q')) 



1 

1 

27n 



e r(0 

d^ 1 



c 



dw 



e " S 







a 
^ dw 
a 



w 
w 



«s(Q') 
«s(Q') 



ln«- 



(3.19) 

(3.20) 
(3.21) 
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which differs from Eq. (3.17) by higher-twist terms, and therefore provides an equally good 
resummation prescription; the difference is in practice very small. 

If we only wish to retain terms which do not vanish as 2; — )■ 1 we may expand 



In - = 1 - z + 0{{1 - zY) 



(3.22) 



with the result 



SBp(^,as(g')) 



1 

2TTi 



1 



c 



dw 



a 



(3.23) 



Equations (3.18), (3.21), and (3.23) provide three alternate definitions of the resummed E func- 



tion which differ by terms suppressed by powers of 1 — 2;. The difference between the first 



two is negligible and the prescription Eq. (3.18) will not be discussed further, but Eqs. (3.21) 



and (3.23) differ by a sizable amount, as we shall see in somewhat greater detail in Sect. 3.2 
Equation (3.23) with C = a is the default Borel prescription of Ref. [131 E 




Figure 6: The damping factor Eq. (3.25) in the Borel inversion of In for various values 
of c= 



When using the Borel prescription, the divergent series Eq. (3.7) is made convergent by 



cutting off its high-order behaviour: indeed, using the Borel prescription Eq. (3.21) to compute 
the inverse Mellin transform of In'' one gets 



27rz 



N+ioo 



dN z-^ In 



N—ioo 



N 



fkiC/a) 



c 



Ck{.z) 

Ini 



(3.24) 
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and similarly if the Borel prescription Eq. (3.23) is adopted instead. Equation 3.24 differs from 
the exact result Eq. (B.4) by the factor fk{C /a) [13j, given by 

7(fc + l,c) 



fk{c) 



(3.25) 



where 7 is the truncated Gamma function. This factor, plotted in Fig. |6j effectively truncates 
the series before the divergence sets in. 



3.2 Subleading terms 



The Borel and minimal prescriptions summarized in Sect. 3J^ differ in the way the high-order 
behaviour of the divergent series is handled. This, as discussed in Refs. [I3l HI], makes in 

practice a small difference unless the hadronic r is close to the Landau pole of the strong 

2 

coupling, tl = 1 — ^-^^j", which is seldom the case, and never for applications at collider 

energies that we are mostly interested in. This is a consequence of the fact that for values of as 
in the perturbative region it is only at very high orders that the effect of the various prescriptions 
kicks in. For example, with the Borel prescription and the "minimal" choice C = 2, if ag = 0.11 
then c ~ 15, so from Fig. [6] one sees that the perturbative expansion is truncated between the 
tenth and twentieth order. 

However, prescriptions may also differ in the subleading terms which are introduced when 
performing the resummation. To understand this, consider the result one gets applying the 
various prescriptions to any finite truncation of the expansion of E(L,as((5^)) powers of 
as(Q^). The minimal prescription then just gives the exact Mellin inverse Eq. (B.6). Because 
this result depends on 1 — z through In 1, in 2; space it generates a series of power suppressed 



terms Eq. (3.22) 



If the Borel prescription is defined according to Eq. (3.21), then the result Eq. (3.24) is 
obtained, which for k small enough that ~ 1 coincides with the minimal prescription, including 
subleading terms. However, with the Borel prescription the z dependence is under analytic 
control: it is entirely contained in the factor ln^~^ - in Eq. (3.21), and it can thus be modified 



at will. Indeed, while the original Borel prescription Eq. (3.21) treats subleading terms as the 



minimal prescription, and thus it coincides with it for finite not too high order truncations of the 



perturbative expansion, using our default Borel prescription Eq. (3.23) instead of Eq. (3.21) only 
leading power contributions in 1 — z are retained. However, this implies that power-suppressed 



contributions in 1/A^ are introduced: indeed, the exact Mellin transforms of Eq. (3.21) and 



Eq. (3.23) are respectively given by 



dzz''-'i:ciz,as{Q'')) 



2TTi 



J 
d^ 



[N^^ - 1] / 
T{N) 



c 



dw 



a 



w 



,«s(Q' 



e [r{N + o r(i + 0J 







dw 
a 



w 



asm 



(3.26) 



which differ by terms suppressed by powers of 
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Hence, the BP Eq. (3.23) on the one hand, and the MP (and the BP Eq. (3.18)) on the other 



hand correspond to two opposite extreme choices in the treatment of subleading terms: in the 
MP, all 1/A^ powers suppressed terms in A'"-space are set to zero, but this leads to 1 — z power 
suppressed terms in z space, while in the BP the opposite is true. The difference between the 



MP and the default BP Eq. (3.23) could thus be taken as a maximal estimate of the impact of 
subleading terms. 

With the Borel prescription it is also possible to construct intermediate, possibly optimized, 
choices of subleading terms. Indeed, in Sect. 2.3| we have noti ced th at in the case of the NLO 



coefficient function Ci(iV), the Mellin transform C}°^(A^) Eq. (2.27) of the logarithmically en- 
hanced 2;-space terms can be brought in better agreement with the full result by inclusion in it 
of some terms suppressed by powers oi z — see Cf^ (A^), defined in Eq. (|2^. This increased 



agreement c an be understood by inspection of the z-space expression of the improved result, 
C}°^ [z] Eq. (2.30): the extra power-suppressed terms which are introduced in it turn out to be 



present in the full result. These terms are in fact of kinematical origin: soft resummation follows 
from the kinematic fact that as 2; — )■ 1 the dependence of partonic cross-sections on z is always 
in the combination Q^(l — -z)^ j23], essentially because this is the upper limit of the integral over 
the energy of radiated gluons for the Drell-Yan process. However, the upper integration limit is 
really 



Q\\-zY 

4:Z 



(3.27) 



so that in fact the resummation produces logarithmic terms of the form In 

This suggests to define an all-order generalization of the approximation (2), by simply 
letting 



ln(l — 2;) — )■ In 



in all resummed expressions. In the Borel prescription, this is easily done: 



1 

27rz 



(3.28) 



(3.29) 



With this choice, the kinematic correction Eq. (3.28) is automatically included to all orders. 
However, beyond 0{as) other subleading terms of the same form but not of kinematic origin 
will in general be present. This choice also arises in a natural way in the context of soft-collinear 
effective theories, and it was adopted in Ref. [7|. 



It turns out that the modified Borel prescription (BP') Eq. (3.29) is closer to the MP than 



the default one Eq. (3.23) because 



In 



1 - z 



so that 



In'^ In i 

Z_ 

Ini 



1 - z 



(l + 0((l-z)^)) 



In'^ ^^ (1 + 0((1 - z)')) 



(3.30) 



(3.31) 
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Equation (3.31) shows that, amusingly, up to terms suppressed by two powers of 1 — z, the 



minimal prescription effectively also performs the kinematic subleading replacement Eq. (3.28) 



though at the cost of also introducing an overall factor ^Jz which is absent in the known pertur- 
bative contributions. 
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Figure 7: The next-to-leading log (NLL) resummed Drell-Yan matched to the fixed next-to- 
leading order (NLO, top row) or next-to-next-to-leading order (NNLO bottom row), with various 
resummation prescription. The right plots show the very large z region. The three alternate 
versions Eq. (3.21), Eq. (3.23) and Eq. (3.29) of the Borel prescription are respectively denote 
BP with log 1/z, BP with (1 - z) and BP with (1 - z)l^. 



In order to assess the impact of these different prescriptions, in Fig.[7]we compare the matched 
result Eq. (3.4) for the Drell-Yan coefficient function in the quark-ant iquark channel obtained 
by including terms up to NLL in the resummed expression, and either up to 0[ol^ (NLO) or 
0(ag) (NNLO) in the fixed-order expansion, with various resummation prescriptions. First, we 
note that the minimal prescription and the Borel prescription Eq. (3.21) (denoted as BP with 
log 1/z in the figure) are essentially indistinguishable (for values of z less than about 0.9, where 
the oscillatory behaviour of the minimal prescription sets in). This is what one expects, since 
they contain the same subleading terms and they only differ in the treatment of the high-order 
divergence. 

However, the Borel prescription Eq. (3.23) is seen to differ by a non-negligible amount from 
the minimal prescription: at large z > 0.3 where the resummation kicks in it is small in compar- 
ison to the size of the resummation itself, while at small z < 0.03, where the resummation just 
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leads to unreliable subleading contributions, it is smaller than the typical higher order correction, 
as it is seen by comparing results matched to the NLO and NNLO. However, in the intermediate 
z region the subleading terms introduced by this prescription are uncomfortably large. 

The improved BP' Borel prescription Eq. (3.29) (denoted as BP with (1 — z)/a/^ in the 
figure) as expected diff'ers less from the MP both at large and at small z] also, it does not 
introduce unnaturally large subleading terms for any value of z. The difference between MP 
and this last version of the Borel prescription is not negligible but small in comparison to the 
size of resummation effects in the region where the resummation is relevant, and it is smaller 
than typical higher-order terms in the region in which the resummation is not relevant. It can 
thus be taken as a reliable estimate of the ambiguity in the resummation. It is interesting to 
observe that the subleading terms which are introduced by the replacement Eq. (3.28) grow at 
small z, and in fact are thus reproducing part of the small-z growth of perturbative coefficient 
functions. The fact that inclusion of these terms is important in keeping the ambiguities of large 
z resummation under control suggests that the large z and small z resummation regions are not 
well separated, and that there might be an interplay between small- and large-z resummation. 



4 Rapidity distribution 

The resummed expression for rapidity distributions was only derived relatively recently in 
Refs. [25| E], confirming a conjecture of Ref. [24]. In this Section, after introducing rapidity 
distributions at fixed perturbative order, we will briefiy review this result, and also compare it 
to a somewhat different expression later derived using SCET methods in Ref. [7]. We will then 
discuss the numerical implementation of resummed results that will be used for phenomenology 
in Sect. [5l 

The hadronic rapidity Y distribution for a Drell-Yan pair of invariant mass Q produced in 
hadronic collisions at center-of-mass energy ^J~s is given by 

(4.1) 



where 



y = Y--\n^, x\ = ^^, a;° = y^e-^, r = ^. (4.2) 

2 X2 s 

The sum runs over all partons in hadrons 1 and 2. To simplify notations, in the following we 
suppress the explicit dependence of the partonic cross-sections doij on and on the factorization 
and renormalization scales fipyfi^. We define 

so that 

1 da dxi dx2 „i 



r dQ'^dV 



Y.I — I — /axi)/|(x2)a. (—.y) . (4.4) 
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Especially in fixed-target experiments, distributions sometimes are given in terms of the 
Feynman xp variable instead of rapidity. The variable a;^ is defined by 



2qL 



(4.5) 



where is the longitudinal momentum of the Drell-Yan pair, and it is related to the rapidity 
Y and the transverse momentum qt by 



1 y^TMlT^) + XF 

r = - log . 

2 ^xj, + 4r(l + ql) - Xp 



(4.6) 



where g|. = q^/Q'^ ■ At leading order = 0, and at NLO is fixed uniquely in terms of Xi? by 
the kinematics, so up to this order the xp and rapidity distributions are simply proportional: 



da 



dQMY 

though at higher orders they will be different. 



NLO ' ^ ^ dQ'^dxp 



(4.7) 



NLO 



4.1 Fixed-order results 

At ncxt-to-lcading order, the rapidity distribution receives contributions from quark-antiquark 
and quark-gluon subprocesses: 



da 



NLO 



rfcT^LO ^ NLO 



dQ'^dY dQ^dY dQ^dY' 
The result is conveniently expressed in terms of new variables z,u, defined by 



Xi 



+ (1 - z)u 
l-{l-z)u 



X2 



-Y 



1 - {l-z)u 
z + [1 — z)u 



with inverse 



e-'^y - z 



X1X2 



u 



^l-z)(l + e- 



X2 



The partonic threshold region Q"^ ^ s = X1X2S corresponds to 2; ^ 1. 
At order the qq contribution is given by 



idaSl 
T dQ^ dY 

where we have defined 



^ dz ^ 
— / duLqq[z,u) 5{1 — z) + -^CpF[z,u) 
z In L 27r 



Lqq{z, U)^J2 ^qqfqi^l) /|(^2), 



(4.8) 



(4.9) 



(4.10) 



(4.11) 
(4.12) 
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the constants Cqq are suitable combinations of the couphng constants of partons to vector boson, 
and 



F{z,u) = {6{u) + 6{l-u) 



/ln(l 



V 1 - ^ 



+ In 



1 - z 



hiz + 1 — z 



+ 



l + z' 



l-z 

The qg + gq contribution is given by 

^Uz 



+ 



1 -M 



2(1 - z) 



(4.13) 



T dQ^ dY 271 



du [Lqg{z, u) G{z, u) + Lgq{z, u) G{z, 1 — u)~\ 



where 



Lqg{z,u) = ^Cqgf^{xi)f^{x2), (q^g), 



(4.14) 
(4.15) 



and 

G{z, u) = 6{u) 



{z' + {l-zf) In 



(1-.) 



+ ln „ 

Z fij 



2z{l-z) 



+ {z' + {l-z)') (^i^ +2z{l 



z) + {l- zfu. (4.16) 



The inclusive coefficient functions, whose form in the quark-antiquark channel and for /ip = 
has already been discussed in Sect. 2J2_, are found integrating over the hadronic rapidity Y . 
We get 



with 



1^ 
T dQ^ 



dQ^ 

^ dz 



6{1-z) + -^CfF,^,{z) 
zn 



2n' 



^ dz 
z 



qg 



T 



gq 



T 



G'int(^) 



and 



: 2Ci{z) + 2\n 
Gint(^) = {z' + (1 - zf) (^In 



{l-zf 



In 



^4 



?>z - -z' 



(4.17) 

(4.18) 
(4.19) 

(4.20) 
(4.21) 
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where Ci{z) was given in Eq. (2.24) and 



^'jW= /"vE'^'j/ZW/K^). (4-22) 



The NNLO rapidity distribution has been computed in Ref. [I], where its lengthy analytic 
expression can be found. 

4.2 Resummation 

As already mentioned, large threshold logs only appear in the quark-quark channel, while the 
contributions from other channels are suppressed by at least one more power of 1 — z as z — )■ 1 



in comparison to the qq contribution. We therefore consider the qq term of Eq. (|4.4|) 
1 da, 



'qq 



T dQ^dY 



where C{z,y) = Cqg{z,y) is the quark- ant iquark coefficient function, which is independent of 
the flavour of the colliding quarks. 

Threshold resummation of rapidity distributions is based on the observation [251 E] (conjec- 
tured in Ref. [2l]) that at large z the coefficient function C{z,y) factorizes as 

C{z,y) = C{z)6{y)[l + 0{l-z)], (4.24) 

where C{z) is the rapidity-integrated coefficient. This is easily proved by rewriting C{z,y) in 
terms of its Fourier transform with respect to y: 

C{z,M)= / dye'''yC{z,y). (4.25) 



The integration range in Eq. (4.25) is restricted by kinematics to \ny/z <y< — \ny/z. Hence, 
one may expand the exponential e**^^ in powers of j/, 

C{z,M)= [ '^^\yCiz,y)[l + Oiy)] = Ciz)[l + 0{l-z)], (4.26) 

Jin 

where C{z) is the rapidity-integrated coefficient function. Hence, C{z,M) is independent of M 
up to terms which as z — )■ 1 are suppressed by powers of |lnz| = 1 — z + 0{{1 — z)"^), and one 
immediately gets the desired factorized form: 

C{z, y)= — e-^''y C{z, M) = C{z) 6{y) [1 + 0{1 - z)] . (4.27) 



oo 



Up to power-suppressed terms we can thus write 
1 da'^"" 



T dQ^dY 



f'^ /' 5^ E /|(-^) 4^ -5'"-) C'-(^)^ (4,28) 
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Because Eq. (4.28) only depends on the rapidity-integrated coefficient function, threshold re- 



summation is simply performed by using for the latter the resummed expressions discussed in 
Sect. E 

Note that resummation of the xp distributions can be performed in the same way. This is 
because the kinematical bounds on Y, 



r I < loe 



r 



translates into 



4 < (1 - rf, 



(4.29) 



(4.30) 



and a similar relation holds at the partonic level. Hence the the above argument holds for the 
xp distribution as well. 



The dependence on parton distributions in Eq. (4.28) can be easily rewritten in terms of the 



differential parton luminosity Eq. (|4.12|) with fixed values of the momentum fractions Xj. To this 

y,„, = ^ln— , (4.31) 



purpose, perform the change of variables 

T 



Z 



with the inverse 



Xi 



X1X2 



1 Xi 

- m — . 

2 X2 



X2 



(4.32) 

z V z 

The line Fcm = Y , selected by the delta function, is contained in the integration region for all 
values of z in 

re^l^l <z<l. 



(4.33) 



Hence, 

1 da' 



T dQ^dY 



dz 



re^\Y\ Z 

^ dz 



C-^(z) / dY,^ 5iY - Fem) Yl fi 

1 



re- 



C-{z)L,,^z 

2\Y\ Z \ ^ 



(4.34) 



where we have used Eq. (|4.12|). The integration range in Eq. (|4.34|) can be extended down to 



z = T because Lqg{z, 1/2) vanishes for z < re^'^' 

1 da'^' 



T dQW 



^ dz /I 
^-^C-{z)L,Jz/- 



(4.35) 



By inspection of Eq. (4.9), we see that, for u = 1/2, xi,X2 depend on z through the ratio t/z. 



Therefore, Eq. (|4.35|) has the form of a convolution product: this greatly simplify its treatment, 

;ement 

(4.36) 



which is then analogous to that of the resummed integrated cross-section, with the replacement 

1~ 
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The matched resummed expression for the Drell-Yan cross-section can finally be written in 
analogy to Eq. (3.4): 



da' 



da 



NPLO 



da' 



dQ^dY 



dQ^dY dQ^dY 



E 

j=0 



1 



d^ da'^' 
MdQHY 



(4.37) 



as=0 



Before turning to the numerical implementation of the resummed result, we discuss briefiy a 
different way of relating the resummation of rapidity distribution to that of the inclusive cross- 
section which has been more recently presented in Ref. [7|. This is based on writing the rapidity 



distribution in terms of the variables z,u defined in Eqs. ( |4.9 , 4.10) 



1 da. 



with 



T (ig2 dY 



C{z,u) 



^ dz 



z 



du Lqg{z, u) C{z, u) 







(9(lna;i, lnx2) 



C{z,y{z,u)) . 



(4.38) 



(4.39) 



d{\nz^ u) 

It is then observed that at NLO the logarithmically enhanced terms in the partonic threshold 
limit 2 — )■ 1 appear as coefficients of the combination 5{u) + 5(1 — m), as one can check by 



inspection of Eq. (|4.13|). At higher orders, logarithmic terms in general multiply non-trivial 

(4.40) 



functions of u, but Eq. (|4.9|) implies that the u dependence of xi, X2 is of order 1 — z, so 

C (z, y{z, u)) = [5{u) + 5(1 - u)] C{z) [1 + 0{l - z)] . 



C{z, u 
Hence 

where 



(9(lna:i, In X2) 



diynz, u) 

1 da'"'' 



dz 

/ -C'"^{z)[L^^{z,0) + L,,{z,l)] 



T dQ"^ dY 



C'''\z) 



(4.41) 



(4.42) 



as it is easily seen integrating Eq. ( 4.40[ ) with respect to u between and 1 and noting that 



It follows that 



9(lna;i, lnx2) 
(9(ln z, u) 

1 da''''' 



du 
dy 



z^ du ( dxi 8x2 8x2 8x1 
T 8y 



r dQ^ dY 



^ dz 



C'^'iz) 



8z 8u 8z 8u 

Lqg{z,0) + Lqq{z,l) 



1. 



(4.43) 



(4.44) 



Equation (4.44) differs by power suppressed terms from the resummed result previously 



derived Eq. (|4.35|), as it is easy to check explicitly. Indeed, using Eq. (|4.9|) and expanding 

sy 

+o{{i-zr) 



Lqq{z, 0), Lqq{z, 1) and Lqq{z, 1/2) lu powcrs of z about 2; = 1 it is easy to check that 

Lqq{z,Q)+Lqq{z,l) ^ ( ^' 



Lqq I Z, 



(4.45) 
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space 



Figure 8: The analytic structure of the integrand of (4.46) and the minimal prescription path. 



Because the difference is suppressed by two powers of 1 — 2; one expects it to be small, and 



indeed we have checked that (using the Borel prescription) the difference between Eqs. (4.35) 



and (4.44) is negligible, and specifically much smaller than the difference between Borel and 
minimal prescription. 

We note that, because Lgg{z,0) and Lqg{z, 1) are not functions of t/z, the form Eq. (4.44) 



of the resummed result does not have the structure of a convolution and thus would require 



a separate numerical implementation. For the same reason, a comparison between Eqs. (4.35) 



and (4.44) using the minimal prescription cannot be performed, because Eq. (4.44) cannot be 
expressed in terms of the Mellin transform of C^'^^{z). We will disregard the form Eq. (4.44) of 
rapidity distributions henceforth. 

4.3 Numerical implementation 

For the sake of phenomenology, an efficient numerical implementation of resummed results using 
the various prescriptions is necessary. Such an implementation was hitherto not available and it 
will be discussed here. 

The minimal prescription involves the numerical evaluation of the complex integral 

1 rln"^'^^ 1 rc+ioo / 1 \ 

-rm? - « U ""^'^^'^ -2) ^'"^"•"'^ 

where the integration path is usually chosen as in Fig. |8] in order to make the integral absolutely 
convergent. However, parton densities obtained from data analysis are commonly available as 
functions of z in interpolated form through common interfaces such as LHAPDF [30], and the 
numerical evaluation of their Mellin transform does not converge along the path of integration 
(specifically, for Re < 0) and must be defined by analytic continuation. The option of applying 
the MP to the partonic cross-section, and then convoluting the result with the parton luminosity 
in momentum space is not viable, because the MP does not have the structure of a convolution: 
the partonic cross-section does not vanish for z > 1, and it oscillates wildly in the region 2; ~ 1. 
This problem is discussed in Ref. [12], where it is handled by adding and subtracting the results of 
the minimal prescription evaluated with a fake luminosity which allows for analytic integration. 
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Another possibility, adopted for example in Ref. p], is to use parton distributions whose 
Mellin transform can be computed exactly at the initial scale. This, however, greatly restricts 
the choice of parton distributions, and specifically it prevents the use of current state-of-the-art 
PDFs from global fits. It is thus not suitable for precision phenomenology. 

The method adopted here, based on an idea suggested long ago [31], consists of expanding the 
function Lqg{z, 1/2) (or Cqq{z) for the inclusive cross-section) on a basis of polynomials whose 
Mellin transform can be computed analytically. We have chosen Chebyshev polynomials, for 
which efficient algorithms for the computation of the expansion coefficients are available. The 
details of the procedure are illustrated in Appendix [C| The obvious drawback of this procedure is 
that it must be carried on for each value of the scale /Xp and, in the case of rapidity distribution, 
for each value of r and Y. 

We now turn to the discussion of an implementation issue which is specific to the Borel 



prescription (BP), and has to do with the choice of the cutoff C. As discussed in Sect. 3.1.2 the 
minimal choice is C = 2; however when C > 1 the integration path in Eq. (3.23) includes values 
of ^ with Re ^ < —1, for which the convolution integral diverges. As discussed in Ref. [H], the 



integral can be nevertheless defined by analytic continuation, by subtracting and adding back 
from Cqqir/z)/ z its Taylor expansion in z around z = 1: the subtracted integrals converge, and 
the compensating terms can be determined analytically and continued in the desired region. 

Here we propose a different method which is numerically much more efficient. The idea is to 
perform the convolution integral with the luminosity analytically, before the complex ^ integral 



of Eq. (3.29). When Re ^ > we can use the identity 



under the integral, because the two ensuing integrals are separately convergent. We get 



(4.47) 



''-^[[l-zf-Xz-iCq.iJ-^ 



dz 



where we have defined 



1 - z 



1 - z 



g{z,T) - Cqg{r) 



Uq,Q-Cqq{r) 



Z2 



(4.48) 
(4.49) 



The function g{z, r) is more easily approximated b y an expansion on the basis of Chebyshev 
polynomials than Cqg (^). Proceeding as in Appendix 



C.2 



we find 



g{z,r) = J2U^-^y 

p=0 



(4.50) 



where n is the order of the approximation, and the coefficients bp = hp{T,jj^) can be computed 
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by numerical methods. We have 



Ida' 



T dQ^ 



d^ 1 



c 



dw 



a 



2m J e r(0 Jo 

5^ 6pB + P + 1, 1 - 1 1 - r") - C^,{t) 
,p=0 ^ ^ 



where 



B(6,a; 1 - r) 



dz z 



a-l. 



1 



\h-\ 



(4.51) 



(4.52) 



is the incomplete Beta function. The function B(6, a, 1 — r) is singular at 6 = 0. In Eq. (4.51), 
the first argument of the B functions in the integrand vanishes for non positive integer values 
of = 0, —1, . . . , in correspondence of zeros of l/r(,^). Thus, the ^ integrand has only a branch 
cut in —w < ^ < 0. This expression is in fact valid for arbitrarily large values of C, hence it can 
also be used in the case C > 2, which requires multiple subtractions if the method of Ref. 
is used. 

The rapidity distribution is obtained by repeating the above procedure with 



11 



L 



<?<? 



(4.53) 



Since L{z, ^) vanishes for z < re the convolution integral is actually restricted to the region 
7-g2|y| < z <1. This difference amounts to replacing r by re^'^': 



1 dal^^ 



TdQ^dY 2Tii J i V{i) 



di 1 



c 



dw 



e <i S 



a 



n 

E 

p=0 



6,B(? + p+l,l-f;l-re^l'"l) 



(4.54) 
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5 Phenomenology 



We will now turn to a phenomenological assessment of the impact of resummation for Drell-Yan 
invariant mass and rapidity distributions. As discussed in Sect. [T| threshold resummation is not 
currently included in the calculations of fixed-target and collider Tevatron Drell-Yan data which 
are included in global parton fits such as MSTW08 [32], NNPDF2.0 [8|, or CTIO [33]; neither it 
has been so far included in predictions for the LHC, such as those of Ref. [3^ . 

We will therefore consider three cases: pp collisions at a center-of-mass energy of 38.76 GeV, 
which corresponds to the case of the experiment E866/NuSea, taken as representative of Tevatron 
fixed-target experiments; pp collisions at the Tevatron collider energy of ^/s = 1.960 TeV; and 
the LHC case, pp collisions at the intermediate energy of ^/s = 7 TeV and at the design energy 
^/s = 14 TeV. For the Tevatron and LHC configurations, we will consider both charged (iu) 
and neutral {££) Drell-Yan pairs, taking into account in the latter case the interference between 
virtual Z and 7. Lepton masses will always be neglected. We will show results for both the 
invariant mass distribution as a function of r = Q'^/s, and for the doubly-differential distribution 
in invariant mass and rapidity as a function of the rapidity for fixed values of r. Specifically, for 
invariant mass distributions we will show results for the /^-factor defined as 

da f 



where 11 p and are the factorization and renormalization scale, respectively. Since we will 
be considering different experimental configurations, the results for the ii'-factors will be shown 
for fixed value of s, with determined by the value of r, = ts. The Born cross-section 

d(T^° , r 

V, 1, 1), which provides the scale for these plots, is shown in Fig. 9 



for LHC at 7 TeV and 



dQ 

Tevatron at 1.96 TeV. 

Our aim will be to assess the potential impact of inclusion of resummation effects on cross- 
sections and their associate perturbative uncertainty, both on experiments which are already 
used for PDF determinations (and thus, potentially, on the PDF extraction itself), as well as on 
future LHC measurements, both for real W and Z production as well as for high-mass 1 TeV 
Drell-Yan pair (relevant for instance as a background to hypothetical Z' production). For each 
observable, we will show fixed-order predictions at leading, next-to-leading and next-to-next- 
to-leading order, and, correspondingly, leading, next-to-leading and next-to-next-to-leading log 
resummed curves. 

All curves will be computed using a fixed (NLO) set of parton distributions. In a realistic 
situation, parton distributions would be refitted each time at the corresponding perturbative 
order; the effect of the perturbative corrections on the hard cross-section is then partly reab- 
sorbed in the PDFs (with fixed experimental data), and the effect on the Drell-Yan process gets 
tangled with the effect on other processes which are used for PDF fitting. Hence, a comparative 
assessment of size of various perturbative corrections on cross-sections and uncertainties, which 
is our main aim here, can only be done with fixed PDF. It should be born in mind, however, 
that our predictions will only be fully realistic when considering the NLO case. 
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Leading-order DY T distribution. Coliider: pp 
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Leading-order DY T distribution. Coiiider: ppbar 
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Figure 9: The invariant mass distribution of charged and neutral Drell-Yan pairs in pp coUisions 
at ^/s = 7 TeV (left) and in pp collisions at ^/s = 1.96 TeV (right) computed at leading order 
with NNPDF2.0 parton distributions at as(M|) = 0.118. 



In order to assess perturbative uncertainties, we will perform standard variations of factoriza- 
tion and renormalization scales, and furthermore in order to assess the ambiguities related to the 
resummation procedure we will compare results obtained with the minimal and Borel prescrip- 
tions, as discussed in Sect.|3| specifically, for the Borel prescription we will use the modified BP' 



Eq. (3.29), which provides a more moderate but realistic estimate of ambiguities as discussed in 



Sect. 3.2, and take C = 2 (see Sect. 3.1.2). Note that we have checked explicitly that also at the 



hadronic level curves obtained using the Borel prescription (3.21), which contains logarithms of 



-, are indistinguishable from those obtained with the minimal prescription, in agreement with 
our discussion in Sect. 13.21 



Other sources of uncertainty will be discussed briefly in Sect. 5.1, where we will provide 
an overall assessment of uncertainties related to the value of the strong coupling and to the 
parton distributions, and then present and evaluate critically the use of scale variation to assess 
perturbative uncertainties. In the remainder of this Section detailed predictions for the three 
cases of interest will be presented. 

All the numerical results of this Section are obtained using the C++ code ReDY (Resummed 
Drell-Yan) [33] . The computation of fixed-order cross-sections in this code relies on the Vrap code 
by L. Dixon [SB], supplemented by soft-gluon resummation described in the previous Sections. 



5.1 Uncertainties 

Theoretical predictions for the Drell-Yan process are affected by a number of uncertainties, 
related to the treatment of both the strong and electroweak interactions. Of course, in a realistic 
experimental situation further uncertainties arise because of the need to introduce kinematic 
cuts, which in turns requires comparing to fully exclusive calculations [37] . Here, we will make 
no attempt to estimate the latter, nor electroweak uncertainties and their interplay with strong 
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corrections (see e.g. Ref. [34j for a recent discussion), and we will concentrate on uncertainties 
related to the treatment of the strong interactions. Before turning to an assessment of the way 
higher order corrections can be estimated from scale variation, we discuss uncertainties related 
to the value of the strong coupling and to the choice of parton distributions (PDFs). 



5.2 Uncertainties due to the parton distributions and as 

The uncertainty due to PDFs is usually dominant in hadron collider processes. Tevatron Drell- 
Yan data are used for PDF determination, so PDF uncertainties here reflect essentially the 
current theoretical uncertainty in knowledge of this process, as well as possible tension between 
Drell-Yan data and other data which go into global PDF fit (which however seems [8] to be very 
moderate). Predictions for the LHC are affected by sizable PDF uncertainties because of the 
need to extrapolate to a new kinematical region, and also, in the case of Drell-Yan, because at 
the LHC, unlike at the Tevatron, one of the two PDFs which enter the leading-order process is 
sea-like. 

PDF uncertainties for the invariant mass distribution of neutral Drell-Yan pairs at ^/s = 



7 TeV are shown in Fig. [lO|as a function of r = Q'^/s. We use NNPDF2.0 PDFs with as(M|) = 
0.118; other PDF sets are expected to give similar results [38] • Because we are using a fixed PDF 
set, the uncertainty does not depend significantly on the perturbative order, or the inclusion of 
resummation. It ranges between 5% and 15% for r < 0.1. For larger values of r the cross-section 
becomes essentially undetermined, because there are no data in PDF global fits to constrain 
PDFs in that region: the few available large-x data are at much lower scale, and the uncertainty 
due to lack of information at very large x > 0.5 contaminates PDFs down to x > 0.1 when 
evolving up to the LHC scale. Note however that the Drell-Yan cross-section at large r > 0.1 
rapidly drops to unmeasurably small values (see Fig. |9]). The fact that PDF uncertainties blow 
up for r > 0.1 implies that data in this region would allow a determination of PDFs in a region 
where they are currently almost unknown; conversely, any signal of new physics in this region 
would have to be validated by measurements in an independent channel (such as for example 
jet production) which provides an independent constraint on the relevant PDFs. 



In Fig. 11 we show the PDF uncertainties for the rapidity distribution of neutral Drell-Yan 
pairs with Q = 1 TeV at ^ = 7 TeV, again using the NNPDF2.0 set with as(M|) = 0.118. 
As in the previous case, the PDF uncertainty does not depend significantly on the perturbative 
order, and it is typically larger than 5%. 

We now turn to the uncertainty due to the value of ag. The current PDG |39] value for as(M|) 
is taken from Ref. [IQJ and it is 0.1184 ± 0.0007. However, this uncertainty seems quite small, 
especially when taking into account the fluctuation in central values between the determinations 
that go into it, and the dependence on the perturbative order of some of them: indeed, current 
recommendations for precision LHC studies from the PDF4LHC group [38] advocate a rather 
more conservative uncertainty estimate. We thus take 

as(M|) = 0.118 ±0.002 (5.2) 

reasonable current range. 
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Figure 10: Invariant mass distribution of neutral Drell-Yan pairs in pp collisions at ^/s = 7 TeV. 
The band corresponds to the PDF uncertainty, using the NNPDF2.0 set with as{M^) = 0.118. 



The impact of this uncertainty on on the Drell-Yan cross-section at ^/s = 7 TeV can be 
estimated from Fig. 12, where we show the effect on the inclusive Drell-Yan cross-section due to 
variation in the range Eq. (5.2) of the value of as{M^) used both in the computation of the hard 
matrix element and the PDF evolution. Results are only shown for r < 0.1, because for larger 
value the PDF uncertainty blows up and results loose significance, as discussed above. Note that 
this full dependence of the physical cross-section on as is in general somewhat different from the 
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Figure 11: Rapidity distribution of neutral Drell-Yan pairs with invariant mass Q = 1 TeV in pp 
collisions at ^/s = 7 TeV. Unresummed results are shown for negative Y and resummed results 
for positive Y. The band corresponds to the PDF uncertainty, obtained by the NNPDF2.0 set 
with as(M|) = 0.118. 



dependence of the hard matrix element alone, because of the dependence on ag of the relevant 
parton luminosity. This total dependence might be larger or smaller according to whether the 
luminosity is correlated or anticorrelated to the value of as, either of which might be the case for 
a quark luminosity, according to the kinematic region [8]. A priori, the size of the uncertainty 
due to variation of as in the matrix element and that due to the dependence on as of the PDFs 
are likely to be comparable: after all, the Drell-Yan rapidity distribution plays a significant role 
in the determination of the PDFs themselves. 

It appears from Fig. 12 that the as uncertainty increases with the perturbative order, but 
it is of similar size at the resummed and unresummed level; at NNLO it is of order of ~ 1.5% 
at LHC energies; we have checked that it is about a factor two larger at Tevatron fixed-target 
experiments. The uncertainty due to ag on rapidity distributions is clearly of comparable size. 

Noting that within the approximation of linear error propagation the PDF and as uncertain- 
ties should be combined in quadrature [H], we conclude that PDF uncertainties are somewhat 
larger than as uncertainties and the combined effect of PDF and as uncertainties is likely to 
be smaller than about 10% but not much smaller, at least in the region in which PDFs are 
constrained by presently available data. Once PDF uncertainties will be reduced due to LHC 
data, it should be possible to keep the combined effect of these uncertainties at the level of 
few percent. Therefore, perturbative accuracies at the percent level are relevant for precision 
phenomenology. 
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Figure 12: Invariant mass distribution of neutral Drell-Yan pairs in pp collisions at ^/s = 7 TeV. 
The uncertainty bands corresponds to a variation of as(M|) in the range 0.116 to 0.120 in the 
hard matrix element and in the parton distributions; NNPDF2.0 parton distributions are used. 



5.3 Perturbative uncertainties: scale variations 

A standard way of estimating unknown higher order perturbative corrections is to vary factoriza- 
tion and renormalization scales. We perform this variation by writing the factorized expression 
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Figure 13: Scale variation grid; the cross-section is computed in correspondence of the purple 
dots. 



Eq. (2.7) as 



aiN,Q')=CiN,i,l)a{N,aM 



^4 



(5.3) 



which is independent of yUp and yu| at the order at which the partonic coefficient a \^N, as{fiD, ^ 

is computed. The residual scale dependence is therefore of the first neglected order, and can be 
used as an estimate of the higher order terms in the perturbative expansion. We vary the two 
scales in the range 





< In 2, 




< In 2, 


In^ 











< In 2, 



(5.4) 



depicted in Fig. [13} which guarantees that both higher-order corrections to the partonic cross- 
section and to perturbative QCD evolution are generated, with the last condition ensuring that 
no artificially large scale ratios are introduced. 

In the sequel, we will perform scale variation of both unresummed and resummed cross- 
sections. The interpretation of results deserves a comment. When performing scale variation 
of a result determined at fixed C(ag), one generates terms of C(ag+^): consequently, the scale 
uncertainty is reduced as one increases the perturbative order. However, terms generated by 
scale variation are proportional to those which are present at the given order: therefore, scale 
variation underestimates the size of higher order corrections when these are enhanced by higher 
logarithmic powers. For instance, scale variation of the 0{as) Drell-Yan coefficient function 
Ci{N) Eq. (2.25) produces terms which at large grow at most as In^ A^, whereas the actual 
0{al) C2{N) coefficient function at large grows as In^A^. Hence, if is so large that these 
terms dominate the coefficient functions and must be resummed to all orders their impact might 
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be rather larger than the scale variation of the fixed-order result may suggest. Nevertheless, 
in this case the scale dependence of the resummed result will still be smaller than that of the 
fixed-order result because the resummed result includes the dominant contributions to the cross- 
section to all orders. 



However, in Sect. 2^ we have seen that there is an intermediate kinematic region in which 
logarithmically enhanced contributions may provide a sizable fraction of the coefficient function 
even though ctg In^ <^ 1: in this case, the resummation improves the fixed-order result in that 
it includes a sizable fraction of the higher order correction, but it still behaves in a perturbative 
way, i.e. terms of higher order in as included through the resummation give an increasingly small 
contribution. If so, the scale dependence of the resummed result may well be comparable to or 
even larger than that of the fixed-order result, because the resummation amounts to the inclusion 
of large terms in the next few higher orders, which are not necessarily seen when performing the 
scale variation of the lower orders. Furthermore, resummation only affects the quark channel, 
while fixed-order scale variation mixes the quark and gluon channels: in an intermediate region, 
the logarithmic terms in the quark channel may be sizable, but with the gluon channel not 
being entirely negligible. In such case, the scale variation is dominated by subleading terms and 
thus we expect the residual scale dependence of the resummed result to differ according to the 
resummation prescription. We will see that this is indeed the case for resummation of Tevatron 
rapidity distributions, with scale variation of resummed results different according to whether 
the Borel or minimal prescription is used. 

5.4 Tevatron at fixed target: NuSea 

We begin by studying the invariant mass distribution of lepton pairs produced by collisions of a 
proton beam of energy E = 800 GeV on a proton or deuteron target, at rest in the laboratory 
(a/s = 38.76 GeV). This is the experimental configuration of the experiment E866/NuSea. 
We first consider the inclusive invariant mass distribution. Results are shown in Fig. Il4| All 



uncertainties shown here and henceforth are due to scale variation as described in Sect. 15.31 As 
expected, the width of the error bands decreases with increasing perturbative order. Note that 
for sufficiently small r the uncertainty blows up, due to the fact that for fixed s the small r limit 
corresponds to low scale: for example, at this energy r = 10"'^ corresponds to Q ~ 1.2 GeV, 



and varying the scales as in Eq. (5.4) the values fin, fip ~ 0.6 GeV are reached. 

Turning now to the resummed results, we note that the numerical impact of resummation is 
large for r > 0.1, while for 0.03 < r < 0.1 it is moderate but sill not negligible. Furthermore, 
starting with the NLL level, the scale uncertainty band for resummed results is dramatically 
smaller than in the case of fixed-order results. This is because scale variation of the LL result 
produces NLL terms which beyond the first few orders are not contained in the fixed-order result; 
starting with NLL these terms are already included in the resummed result. It is interesting to 
note that in the case of the resummed cross-section (with both prescriptions) the NNLL band is 
almost entirely contained in the NLL band, while the fixed-order NLO and NNLO error bands 
are only marginally compatible with each other. The ambiguity in the resummation, estimated 
from the difference between Borel and minimal prescription, is not negligible, but smaller than 
the scale uncertainty; moreover, it is more evident at small r, since the different subleading 
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Figure 14: Invariant mass distribution of neutral Drell-Yan pairs in pp collisions at ^/s = 
38.76 GeV. 



terms give a larger contribution in that region. 

The experiment E866/NuSea [1] has measured the distribution in xp Eq. (4.5) of lepton pairs 
with an invariant mass Q = 8 GeV. The E866 data are displayed in Fig. 15, superimposed to the 
QCD prediction and the corresponding scale uncertainty. The distribution is symmetric about 
Y = 0; the curves shown for Y < refer to fixed-order calculations, and those with y > to 
resummed results. The data agree with the NLO calculation because these data were included 
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DY rapidity distribution. Coliider: pp Subprocess: Z+gamma 




Figure 15: Rapidity distribution of neutral Drell-Yan pairs of invariant mass Q = 8 GeV in pp 
collisions at ^/s = 38.76 GeV; E866 data are also shown. 



in the determination of the NLO PDFs that we are using. 

The impact of the resummation is small but not negligible: for instance the difference between 
NNLO and NNLL is about half of the difference between NNLO and NLO. Furthermore, the scale 
uncertainty of resummed results is somewhat smaller than that of the unresummed ones. This is 
consistent with the observation that for this experiment r = 0.04, which, as discussed in Sect. [2| 
is in the region in which resummation is relevant. However, the difference between resummed 
results obtained using the Borel and the minimal prescription is almost as large as the size of the 
resummation itself: in fact the NLL Borel results is somewhat lower than the NNLO one, while 
the NLL minimal prescription result is a bit higher. Hence we conclude that the overall impact 
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of the resummation on these data is essentially negligible. A large and negative NLL resummed 
correction to the NLO result was claimed in Ref. [6j, using the minimal prescription, but we do 
not confirm it: we find a positive and rather smaller correction. The result of Ref. |6] was first 
criticized in Ref. [7]. Our result with the Borel prescription is in good quantitative agreement 
with Ref. [7j; however, the minimal prescription gives a somewhat larger correction, though still 
positive. 



5.5 Tevatron collider 

We now turn to Drell-Yan production in pp collisions at a center-of-mass energy y/s = 1.96 TeV. 
Results for the invariant mass distribution of neutral and charged Drell-Yan pairs in this config- 



uration are shown in Fig. 16 
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Figure 16: Invariant mass distribution of neutral (left) and charged (right) Drell-Yan pairs in pp 
collisions at y/s = 1.960 TeV. 



The behaviour of these curves is similar to that seen in the case of NuSea, Fig. 14, but 
with the impact of the resummation yet a bit smaller, as one would expect both because of the 
higher energy and because of the collider configuration, as discussed in Sect. 2.2 (in particular 
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Fig. [3]). Interestingly, even when the resummation has a very small impact, it still leads to 
a non-negligible reduction of the uncertainty: this is consistent with the expectation that for 
these medium-small values of r resummation is in fact a perturbative correction, as discussed in 
the end of Sect. |2.3 Note that in these plots the smallness of leading-order uncertainty bands 
when r ^ 0.002 (i.e. ^ 100 GeV^) is due to the fact that the scale dependence of the parton 



luminosity, to which the leading-order cross-section is proportional (see Eq. (4.11)), is stationary 
at this scale. 
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Figure 17: Rapidity distribution of neutral Drell-Yan pairs of invariant mass Q = Mz in pp 
collisions at y/s = 1.96 TeV (the contribution of virtual 7 at the Z peak is included). 



We now turn to rapidity distributions, shown in Figs. [17] and [18] for neutral Drell-Yan pairs 
of invariant mass Q = Mz and Q = 200 GeV respectively, and in Fig. [19] for charged Drell-Yan 
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Figure 18: Rapidity distribution of neutral Drell-Yan pairs of invariant mass Q = 200 GeV in 
pp collisions at y/s = 1.96 TeV. 
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Figure 19: Rapidity distribution of charged Drell-Yan pairs of invariant mass Q = Mw in pp 
collisions at y/s = 1.96 TeV. 
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pairs of invariant mass Q = Mw- The impact of the resummation is now very small, as one 
would expect given the smallness of the relevant values of r. However, interestingly, resummed 
uncertainty bands are systematically smaller than the unresummed one, with the resummation 
ambiguity (i.e. the difference between minimal and Borel results) now essentially negligible. 
Hence, even in this small r region the resummation leads to perturbative improvement, while 
behaving of course as a perturbative correction. Again, note that the smallness of leading-order 
uncertainty bands is due to the fact that the scale here is close to the stationary point of the 



scale dependence of the parton luminosity already seen in Fig. 16 



We can compare these uncertainties to that of typical current data, thanks to recent measure- 
ments at the Tevatron. Specifically, the rapidity distribution of e^e~ pairs with invariant mass 
in the range 66 GeV < Q < 116 GeV has been recently measured by the CDF collaboration |32]. 
In principle, the data should be compared with the theoretical prediction for the full process 



p + p + e + X. 



(5.5) 



For values of Q close to the Z mass, however, a good approximation is provided by the Breit- 
Wigner approximation, which amounts to assuming that the amplitude is dominated by Z 
exchange, and takes into account the finite width Tz of the Z boson: 



da{T, Y, Q2) 



2MzT, 



1 da: 



dQ'^dY (Q2 _ M|)2 + M|r| 27r dY 



(5.6) 



where Fe+g- is the Z decay width into a lepton pair, and daz is the differential cross-section for 



the production of a real on-shell Z boson. Eq. (5.6) gives 



rfa(r,F,g2 



M|F| 



c/a(r, y, M| 



dQ^dY 



MIY + M|F| 



and therefore 



where 
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is just a y-independent multiplicative factor. 
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(Q2 _ M|)2 + M|F| 



(5.7) 



(5.^ 



(5.9) 



In Fig. 
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we show the CDF data |12], corresponding to an integrated luminosity of 2.1 fb^^, 
compared to the NNLO QCD prediction with the inclusion of threshold resummation at NNLL, 
using both Borel and minimal prescriptions. The comparison shows an excellent agreement in 
shape between the data and the theoretical curves; there is clearly a mismatch in normalization of 
the total cross-section, which is however consistent with the sizable 6% normalization uncertainty 
on the cross-section. Also in this for the NuSea experiment, this simply reflects the fact 

that these data are used in the determination of the PDFs that we are using, and the difference 
between the NLO expression used in PDF fitting and the NNLO one shown here is much smaller 
than the experimental uncertainties. 
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Figure 20: Rapidity distribution of Z bosons in collisions at a/s = 1.96 TeV (the contribution 
of virtual 7 at the Z peak is included). Data are taken from |42j. The smaller black uncertainty 
bands are statistical only, while the wider green bands also include normalization uncertainties. 



A similar comparison can be performed for the asymmetry, defined as 

^w{Y) = -f^ (5.10) 

daw+ daw- 

dY ^ dY 

also measured by CDF [13]. In this case, normalization uncertainties cancel in the ratio. In 
Fig. 21 we show the measured CDF data compared to the QCD prediction at NNLO and 
resummed NNLO+NNLL (Borel and minimal prescriptions). 

Clearly, the accuracy of present-day data is insufficient to appreciate the effect of NNLO or 
resummation correction, and it is rather comparable to the difference between LO and NLO 
predictions, which can thus be barely appreciated. However, an improvement of statistical 
accuracy by an order of magnitude would be sufficient for NNLO and resummation corrections 
to become significant. The normalization uncertainty has a negligible effect on the shape of the 
distribution and therefore it does not affect this conclusion 



5.6 LHC 

We now consider predictions for Drell-Yan production at the LHC, both at 7 TeV and 14 TeV. In- 



variant mass distributions for both neutral and charged Drell-Yan pairs are shown in Figs. 22 , 23 
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Figure 21: asymmetry in pp collisions at ^/s = 1.96 TeV. Data are taken from 



While the impact of fixed-order perturbative corrections is unsurprisingly similar to that at the 
Tevatron collider shown in Fig. 
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interestingly the reduction in uncertainty obtained thanks 
to the resummation is more marked at the LHC, consistent with the expectation (recall Sect. [2| 
that the effect of the resummation is somewhat more significant at a pp than at a pp collider. 
Moreover, the consistency of the NLO error band with the NNLO prediction is improved by the 
inclusion of resummation. Of course, as in the case of Tevatron, for realistic values of r < 0.1 the 
impact of the resummation is mostly on the uncertainty but very small or negligible on central 
values, so the resummation is behaving as a perturbative correction. 

Turning to rapidity distributions, we present results for the following observables: 



neutral Drell-Yan pairs with invariant mass Q = I TeV, Fig. 24 



neutral Drell-Yan pairs with invariant mass Q = Mz, Fig. 25 



positively charged Drell-Yan pairs with invariant mass Q = Mw, Fig. 26 



negatively charged Drell-Yan pair with invariant mass Q = M\y, Fig. 27 



The first case corresponds to r ~ 0.02 
200 GeV at the Tevatron shown in Fig 



comparable to the case of a final state with Q = 
As in that case, we clearly see an improvement in 
uncertainty (with small resummation ambiguities) when going to the resummed result, though 
also in that case the effect on central value is moderate. On the other hand, the other cases 
correspond to very small values of r and indeed in this case the uncertainty on resummed 



47 




Figure 22: Invariant mass distribution of neutral Drell-Yan pairs in pp collisions at ^/s = 7 TeV. 
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Figure 23: Invariant mass distribution of positively (left) and negatively (right) charged Drell- 
Yan pairs in pp collisions at ^/s = 7 TeV. 
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results is larger than that on unresummed ones, indicating that whatever effect is induced by 
the resummation is related to the inclusion of terms which are not dominant in this region. This 
is also reflected in a sizable difference between Borel and minimal results. 
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Figure 24: Rapidity distribution of neutral Drell-Yan pairs of invariant mass Q = 1 TeV in pp 
collisions at a/s = 7 TeV. 



Finally, in Figs. 28 33 results for this same invariant mass and rapidity distributions for the 
LHC at 14 TeV are shown. The general behaviour is the same as that at the lower energy. 
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Figure 25: Rapidity distribution of neutral Drell-Yan pairs of invariant mass Q = Mz in pp 
collisions at y/s = 7 Te V (the contribution of virtual 7 at the Z peak is included) . 



51 



DY rapidity distribution. Collider; pp Subprocess: W+ 



LO 
NLO [ 
NNLO s 
Borel LL+LO 
>(;el NLL+NLO [ 
Bo^iNLL+NNLO ! 




DY rapidity distribution. Collider; pp Subprocess; W+ 



LO 
NLO [ 
NNLO 1 

X Minimal LL+LO [ 
../linimal NLL+NLO 
LL+NNLO ! 




Figure 26: Rapidity distribution of positively charged Drell-Yan pairs of invariant mass Q = M\y 
in pp collisions at ^/s = 7 TeV. 
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Figure 27: Rapidity distribution of negatively charged Drell-Yan pairs of invariant mass Q = My/ 
in pp collisions at ^/s = 7 TeV. 
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Figure 28: Invariant mass distribution of neutral Drell-Yan pairs in pp collisions at ^/s = 14 TeV. 
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Figure 29: Invariant mass distribution of positively (left) and negatively (right) charged Drell- 
Yan pairs in pp collisions at ^/s = 14 TeV. 
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Figure 30: Rapidity distribution of neutral Drell-Yan pairs of invariant mass Q = Mz in pp 
collisions at y/s = 14 TeV (the contribution of virtual 7 at the Z peak is included). 
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Figure 31: Rapidity distribution of Z bosons in collisions at ^J~s = 14 TeV (the contribution 
of virtual 7 at the Z peak is included). 
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Figure 32: Rapidity distribution of positively charged Drell-Yan pairs of invariant mass Q = My/ 
in pp collisions at ^/s = 14 TeV. 
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Figure 33: Rapidity distribution of negatively charged Drell-Yan pairs of invariant mass Q = My/ 
in pp collisions at ^/s = 14 TeV. 
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6 Conclusions 



We have presented a discussion of the theory and phenomenology of threshold resummation of 
rapidity distributions, with special regard to present-day hadron colliders. On the theoretical 
side, our main result is that we have provided a way to estimate the value of r = Q"^ j s at which 
resummation of threshold logarithms is expected to provide an improvement over fixed-order 
calculations. This result has been accomplished through a determination of the relevant partonic 
center-of-mass energy, whose distance from threshold determines the impact of resummation. 
This estimate relies on the singularity structure of the anomalous dimensions which drive 
evolution of parton distributions in perturbative QCD. Using this technique, we have shown that 
resummation is expected to be relevant down to fairly small values of r ~ 0.003 at pp colliders, 
and r ~ 0.02 at colliders. 

At the phenomenological level, we have shown how to use different versions of the Borel 
resummation prescription and their comparison to the minimal prescription as a means to assess 
the ambiguities related to the resummation, and in particular to the treatment of the subleading 
terms. The application of these tools to the Tevatron and LHC has shown that resummation is 
relevant for the production of states of mass as light as W and Z vector bosons at the Tevatron, 
and for the production of heavy dileptons of mass in the TeV range at the LHC: in all these cases 
threshold resummation leads to a significant improvement in perturbative accuracy. The impact 
of resummation on Tevatron fixed-target rapidity distributions is less clear, in that, despite being 
larger, the effect of resummation may be marred by its ambiguities. 

Our general conclusion is that the impact of threshold resummation at hadron colliders for 
T > 0.01 is comparable to that of NNLO corrections, and it should thus be included both in the 
determination of parton distributions and in precision phenomenology, though care should be 
taken in also estimating carefully the ambiguity which is intrinsic in the resummation procedure. 

Acknowledgements: We thank M. Grazzini and A. Vicini for discussions and S. Catani for 
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A Saddle point 



We collect here some results used in Sect. [2] on the saddle-point approximation to Mellin trans- 
form integrals of typical parton distributions. 

First, we use the saddle-point method to show that the large- behaviour of a Mellin trans- 
form describes the original function of x in the large-x region. To see this, consider the Mellin 
transform of a generic function f{x), 



1 

N-l 



f{N)= / dxx^'-'fix); (A.l) 
Jo 

we assume that f{x) > in the range < a; < 1. For real, f{N) is a decreasing function of 
A^, because the area below the curve x^~^f{x) obviously decreases as A^ increases. 
The inverse Mellin transform is given by 



1 



N+ioo 



f(x) = — dNx-''fiN), (A.2) 

with N larger than the real part of the rightmost singularity of f{N). We estimate the range of 
A" which gives the dominant contribution by saddle-point: we rewrite Eq. (A.2) as 

-| pN+ioo 

f{x) = — ciATe^^^'^) (A.3) 

where 

E(x,N) = N\n- + \nf(N). (A A) 

X 

For appropriate f{N), as in the case of the physical cross-section, E{x, N) has a unique minimum 
on the real A^ axis at A^ = A^o, with 

E'{x,No) = \n- + l^ = 0, (A.5) 

where a prime denotes differentiation with respect to A^. The inversion integral is dominated by 
the region of A^ around No, and can be approximated by 



2 



1 fNo+ioo „ 

fix) ^ — / rfiVe^(^'^°)+'^^(^-^°) 

1 f,E{x,No) 

(A.6) 



V2^y/E"{x,No) 

(after the variable change N = Nq + U and gaussian integration). When x — )■ 1, the slope of the 
straight line A^ In decreases, and the position A"o of the minimum is pushed to larger values. 
We conclude that the behaviour of f{x) for x close to 1 is determined by the large- A^ behaviour 
of f{N), as we set out to prove. 
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As an explicit example, we consider 



/(x) = x"ln^-, 

X 



(A.7) 



with a and (5 positive constants. The function f{x) Eq. (A.7) vanishes both at x = and a; = 1, 
and has a maximum at 

X = e"° (A. 8) 

which approaches 1 as /3/a; — > 0. We find 



r(/3 + i) 



(A.9) 



Let us now compute f{x) using the saddle-point approximation to invert the Mellin transform 
Eq. (IXgl). We have 



E(x, N) = N\n- + In r(/3 + 1) - (/3 + 1) ln(A^ + a) 

X 

E'{x, N)=\n-- ^ ^ ^ 



E"{x,N) 



X N + a 
(3 + 1 



(N + a) 



(A.IO) 
(A.ll) 
(A.12) 



From Eq. (A.ll) we find that the saddle point is at 

f3 + l 



Ini 

X 



a, 



(A.13) 



which is, as expected, an increasing function of x. Furthermore 



X 



E"{x,No 



Therefore, Eq. (A. 6) gives 



/3 + 1 X 



fix) ^ ir(/3)x"ln^ 



X 



where 



Ki(3) 



r(/3 + 2) 



(A.14) 
(A.15) 

(A.16) 
(A.17) 



VWTl) exp [{/3 + 1) ln(/3 + 1) - (/3 + 1)] 

We see that the saddle-point approximation in this case gives the exact result, Eq. (A.7), up to 
a normalization factor, which however approaches 1 for large values of /3 (the denominator of 
Eq. (A.17) is just r(/3 + 2) in the Stirling approximation). Indeed, one can check that the Taylor 
expansion of E{x, N) to second order around Nq is increasingly accurate as /3 — ?■ oo. 
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It is important to observe that the argument which shows that a Melhn transform f{N) 



Eq. (A.l) is a decreasing function of does not apply if f{x) in Eq. (A. 2) is a distribution 



rather than an ordinary function. An obvious counter-example is the Dirac delta function, 
which has an A^-independent Mellin transform: 



dxx^-^6{l-x) = 1. 



(A.18) 



It is not difficult to find examples of distributions whose Mellin transforms even increase as A^ 
increases. Let us consider for instance 



d{x) = [D{x)]^, 

where D{x) is a positive function, at most as singular as (1 
definition of the plus prescription, we have 



x)' 



m x 



d{N) 



dx (x 



N-l 



1)D{ 



x) 



X 



N-V 



(A.19) 

= 1. With the usual 
(A.20) 

D(x) increases as A^ 



In 



k 1 

1-3 

1-x 



which is finite under the above assumptions. The area below ( 
increases, and therefore |rf(A^)| grows with A^. 

In fact (see also Appendix B) this is the typical situation for the distributions 

which characterize the large-x and thus large- A^ behaviour of soft-gluon emission cross-sections 
after the cancellation of infrared poles, and whose Mellin transforms behave as In'^^^ A^ when 
N ^ oo. The coefficient of this logarithmic rise need not be positive, given that these are 
partonic cross-sections obtained after subtraction of infrared poles and factorization of collinear 
singularity. 

Of course, the physical cross-section must be an ordinary function, so in the factorized 



expression Eq. (2.1) the decrease of the parton luminosity at large A^ always offsets the increase 
of the partonic cross-section a. Indeed, because the partonic cross-section rises at most as a 
power of In A^ as A^ — )■ oo, it is easy to show that a sufficient condition for the Mellin inversion 



integral Eq. (2.5) of the factorized cross-section Eq. (2.1) to exist is that the parton luminosity 



C{z) vanishes at large z at least as a positive power of {1 — z). 



B Mellin transformation 



We collect here some results on the calculation of Mellin transforms of functions and distributions 
which typically appear in the perturbative coefficients. We start from the identity 



dz z 



ln« 



/ dz {z^-^ - 1) ln«-^ - = r(0 (A^^^ - 1) (B.l) 



(the plus prescription is not necessary if Re^ > 0, but its presence makes the l.h.s. well defined 
even when ^ = 0). It follows that 

-1 



A^-« = 1 + 



1 

m Jo 



dz z 



N-l 



ln« 



-il 

z 



(B.2) 
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and therefore 



1 
27ri 



N+ioo 



dNz-^N-^ = S{l-z) + 



N—ioo 



In' 



;-l 1 



no 



(B.3) 



with > 0. Taking k derivatives with respect to ^ at ^ = we get 

"N+ioo -| 



1 

27ri 



dN z-^ In 



N—ioo 



N 



6ko6{l - z) + 



d' ln«-i I 




de m 





(B.4) 



Note that the result Eq. ( |B.4[ ) is a distribution, consistently with the fact that In'^ is an 
increasing function of N . In other words, the inversion integral does not exist in the ordinary 
sense, but only in the sense of distributions. 



An equivalent form of Eq. (B.4) is 

»7V+ioo 



1 

2tH 



dN z-^ In' 



N—ioo 



N 



4o5(l -z) + 



2Txi 



di ln^-^1 



^.+1 r(e) 



(B.5) 



where the integration contour is any closed curve in the complex plane ^ with the origin ^ = 
inside. This form proves useful in the formulation of the Borel prescription for the Mellin 
inversion of divergent series in In A^. 

A second equivalent form of Eq. (B.4) is obtained by computing the derivative in the r.h.s. ex- 
plicitly: 



1 
2Tii 



N+ioo -1 

dNz-"" ln*=- 



N—ioo 



where 



A(0 



4o^(l - z) + 



1 



In 



In In 



z j=l V/ 



1 



(B.6) 



(B.7) 



We see that the terms in the sum over j are decreasingly important as z — ?■ 1. 

Threshold logarithms typically appear in perturbative coefficients through the Mellin trans- 
form 



h{N) 



f dzz""-^ 


'\n\l-zy 


Jo 


1-z 



We have 



where 



Hence 



d'' 



I^(N)=—^F{N,0 



F{N,0 



dz (z 



N-l 



1){1-Z) 



r(Ar)r(i + o 
r(Ar + o 



Ik{N) 



k + 



k+1 
j=0 



(B.8) 
(B.9) 
(B.IO) 
(B.ll) 
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C Chebyshev polynomials 



In this Appendix we recall the definition and the main properties of Chebyshev polynomials. 
The Chebyshev polynomials 

i 

T,{z) = Y,T^kz\ (C.l) 

fc=0 

are defined in the range z G [—1,1] recursively by 

To(^) = 1 (C.2a) 

Ti{z) = z (C.2b) 

T2{z) = 2z^ -I (C.2c) 

Ti{z)=2zTi_^{z)-Ti^2. (C.2d) 

A generic function G{u) can be approximated in the range u G [Mmin5Wmax] by its expansion 



on the basis of Chebyshev polynomials (C.2), truncated at some order n: 



G{u) ^ -| + 5Z Q Ti{Au + 5), (C.3) 

where 

^max ^min ^max ^min 

Simple numerical algorithms for the calculation of the coefficients Cj are available (we have used 
the routines of the gsl). 



Simple algebra leads to 



G{u)c^^Ck{Au + Bf (C.5) 



where 



/c=0 



Cfc = -y^fcO + ^Q^jfc- (C.6) 

i=k 



C.l Minimal Prescription 

We have shown in Sect. |3.1| that the minimal prescription can be conveniently implemented by 
means of an analytic expression for the Mellin transform of the luminosity C{z) (which can be 
either Cqq{z) in the case of invariant mass distributions or Lqq{z, 1/2) in the case of rapidity 
distributions). This can be obtained by expanding C{z) on the basis of Chebyshev polynomials, 
truncated at some finite order n, and then taking its analytical Mellin transform. The luminosity 
itself, however, is very badly behaved in the range (0, 1): it is singular at 2; = 0, and varies by 
orders of magnitude in the range 

fl for the rapidity- integrated cross-section 

I I for the rapidity distribution 
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It is therefore convenient to expand a regularized function 

F{z) = j-^^^C{z). (C.8) 

Values of (3 and 5 in the range 3 -i- 7 are normally suited to make F{z) smooth enough to be 
approximated by a reasonably small number of Chebyshev polynomials. Equation (C.4) gives 

A=—, B = -l. (C.9) 

and the approximation is 

n y \ k n 

F{z) = j2ckh—-i] =J2^p''^ (c.io) 

k=0 ^ ' p=0 

where 



2 Y^fk 



(c.u) 



The luminosity is easily recovered: 



p=0 j=0 ^-^^ p=0 



where the last equality holds for S integer. It is now immediate to obtain the Mellin transform: 

p=0 j=0 

Alternatively, one may introduce the variable 

Z = ^n.axe" (C.14) 

and appoximate the function 

F{u) = ZC{Z) = Zinax e"£ (2:maxe") (C.15) 

by an expansion on the basis of Chebyshev polynomials. The variable u ranges from — oo 
to when < z < Zj^^ix] however, for practical purposes one only needs the luminosity for 
z > Zmm = -ZmaxC"™". Wc havc thcrcforc 



and 



A = , 5 = 1, (C.16) 

n , V A; 

F(u) = Y,cA\-2—\ . (C.17) 

i._n V "^min / 



fc=0 
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We can now reconstruct C(z) through the replacement u = In We get 

p=0 k=p ^ ^ 

The Melhn transform is computed using the result 



'max 







We obtain 



where 



(C.l^ 



Zmax ^^'^ (A^ - 1)P+1 ■ {y-'-i-J) 



dzz'^-'Ciz) = zl-; g (C.20) 

In practice, we have found that the second method is to be preferred for small values of r, 
r < 0.1, while the previous one works better for r > 0.1. 



C.2 Borel Prescription 



In this case, we look for an approximation of the function (7(2;, r), Eq. (4.49), as a function of 
z G [r, 1] . We have 



g{z, r) = J2ck {Az + Bf = Y, K (1 - 

k=0 p=0 



where 



and 



k=p ■' 



Ck 



A 



1-t' 

Note that A + i? = 1 in this case. Therefore 



B 



1 + r 
1-t' 



(C.22) 



(C.23) 



(C.24) 



-2 fk 



1 - r 



k=p 



P 



Ck- 



(C.25) 



In the case of the rapidity distributions, the variable z is in the range z G [re^' ^ ' , 1] ; therefore 



l-re^\Y\ 



2 Vj^/'k 



k=p 



P 



Ck- 



(C.26) 
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D Resummed cross-section for the Drell-Yan process 

In this Appendix we give the exphcit expressions of the functions Qi which appear in the re- 
summed Drell-Yan cross-section, Eq. (3.2). We have 



9iW 



1 + asgol + algo2 + 0{a 
2Ai 



PI 



[(1 + A) ln(l + A) - A] 



(D.la) 
(D.lb) 



[A-ln(l + A)] + 



/3o 



ln(l + A)(ln^-27, 



- Aln 



,21 



+ 



ln2(l + A) + ln(l + A) - A 



^3(A) 



PI 



(D.lc) 



+ A 



+ 



+ 



Ai(5l ln(l + A) 
2(31 1 + A 



l + -ln(l + A) 



+ 



ln(l + A) 



1 + A 



A, 



1 + A 
D2 



A 

iTa 



ln(l + A) 



Ai , A2 



A 



A2 , , fJ'l Ai 
2(3i ^ 4/3o 



A In- 



1 + A 

2 /fp _ 



+ 



ln(l + A) 



2Pi 
A 



1 + A 



In^ 



1 + A 



Q 



In 



0! 



21 



with 



d 



din fx' 



-asiiJ) = -(3oal (1 + (3ias + (32al + . ■ ■) . 



(D.ld) 



(D.2) 



The coefficients gok can be found in |15], but without scale-dependent terms. Their full expression 
is given by 



9oi 



902 



Cf 
167r2 



4C2 - 4 + 27^ + 
511 



-27, 



(D.3a) 



+ C^ 
+ n/ 
+ In' 



4 



552 

198 C2 - 60 C3 + ^ C2 - 128 7' + 128 7^(2 + 32 7 
5 



Q_ 



1535 376 604 
^2~ ^ IT ^ ~9~ ~ 

127 64 8 224 

- y C2 + g Cs - ^ 7^ 

2 r 



92 



1616 
~27 



■1e 



536 2 2 > 

56 7eC3 + ^ 7b - 16 7bC2 + — 7 



176 



Ci. (327^ - 487^ + 18) + Ca 

68 



80 2 32 3 

Y7b- Y7b 

44 



7^-11+72/ 2 



;7b 



+ In 



Cf (48C3 + 72C2 - 93 - 1287^(2 + 1287^ + 487I - 647I) 

536 88 



+ C^(^-24C3 



-C2 + 167i5C2 



9 



-1e 



-1e 



TT 



+ ^/ ( y C2 



4C2-4 + 272 + 



34 80 



16 2 



27^ In 



In 



The 

/32 



coefficients appearing in the previous functions are 

_ IICa - 4:Tf _ 33 - 2nf 

~ mr ~ 127r 

_ 17^2 - (IOCa + 6Cf)Tf Uf _ 153 - 19nj 



247r2/3o 



(33 - 2n/)27r 



1 



2857 3 / 2 205 1415 ^ 



^2 
^3 



(33 - 2n/)327r2 

Cf_± 
7r 37r 



2857 



5033 



325 



Do 



Cf_ 
27r2 

47]-3 
+ 

Cf 
167r2 



^2 



67 




18 




245 
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10 



9 



201 - lOn 



/ 



C2 + ^ Ca 
6 



7 



11 



277r2 



+ - 



1 
3 



(D.3b) 

(D.4) 
(D.5) 



Tin) 



(D.6) 
(D.7) 
(D.8) 



55 
24 



2Cz]CFnf 



1 



209 10 ^ . ^ , ^ 



Ca 



1616 88 



27 



9 



■vr^ + 56C3 + 



224 



16 



(D.9) 
(D.IO) 



D.l Matching 

Here we compute the terms to be subtracted in the resummed result in order to avoid double 



counting. Start from Eqs. (|D.l|) and expand them in powers of their argument A: 

1 



^?i(A) 
^72 (A) 
^73(A) 



/3o 



A' - ^X' 



A, 



ln^-27,')A+l — 2 + 
l4 J v2/3o 2/3o 



27^. - In 



A=^ + 0(A^ 



(D.lla) 
(D.llb) 



A. 2 



/3o 



(7^ + C2) 



r7i5 + ^t;^ I A 



PI 
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+ 



—7s In — + —2 In — - — I In 



ln^|j)A + 0(A^). 



(D.llc) 



The Taylor expansion of g'o(Q;s) exp5(A,Q;) is 

5o(as) exp5(A,a) = (1 + as^oi + «s5o2 + • • .)e"s'5i+ai52+... 

= 1 + (5i + ^oi)as + Y + ^2 + 5i ^01 + 502^ «s + 0(«s)- 

where (note that, since A = a In it can be seen as an expansion in A with X/a fixed) 



(D.12) 



oisSi = 



alS2 



A, /A g2 

— - + In — - 27, 



A 



^1 A 

3Ai^ ■ \2Pl ■ 2/3o 



+ 



^2 ^A, ( Q2 
+ — I 27b - In — 



^1 



+ — 7Bln 
^0 /^l 



+ 



A, 



g2 



In-^ - 



4/5o 



(D.13) 



+ ^ -^(7. + C2)-^7. + ^ 



ln2 ^ - 



(D.14) 
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